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ABSTRACT 


The  equilibrium  and  stability  of  a  gas-vapor  bubble 
in  a  liquid-gas  solution  at  negative  pressures  is  investigated. 
The  bubble  is  considered  to  exist  at  the  center  of  rotation 
of  a  thin  tube  filled  with  the  solution.  The  centrifugal 
field  generated  by  the  rotation  of  the  system  produces  a 
state  of  tension  in  the  solution  surrounding  the  bubble. 

The  equilibrium  shape  of  the  bubble  is  shown  to  be 
approximated  by  a  spheroid  for  which  the  plane  section  is  an 
ellipse.  However,  for  negative  pressures  in  the  solution  the 
bubble  is  very  close  to  being  spherical.  The  stability 
analysis  shows  that  the  equilibrium  state  of  the  bubble  is 
unstable,  and  the  equilibrium  size  represents  a  critical 
value.  A  bubble  less  than  the  critical  size  will  spontan¬ 
eously  collapse,  while  a  bubble  greater  than  the  critical 
size  will  spontaneously  grow. 

Consideration  is  then  restricted  to  a  dilute,  liquid- 
gas  solution.  The  critical  radius  is  expressed  as  a  function 
of  the  properties  of  the  liquid  phase  outside  the  center  of 
the  bubble.  The  dependence  of  the  critical  radius  on  the 
pressure,  concentration,  and  temperature  is  investigated  for 
the  specific  case  of  a  dilute  solution  of  air  in  water.  The 
critical  radius  is  found  to  decrease  as  the  magnitude  of  the 
negative  pressure  increases.  The  presence  of  the  dissolved 
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gas  results  in  a  small  reduction  in  the  critical  radius  at 
negative  pressures. 


v 


ACKNOWLEDGEMENT 


I  wish  to  thank  my  supervisor,  Dr.  T.  W.  Forest, 
for  the  guidance  and  assistance  he  provided  throughout  the 
course  of  this  work. 


vi 


TABLE  OF  CONTENTS 


CHAPTER  PAGE 

1  INTRODUCTION  .  1 

2  THE  EQUILIBRIUM  STATE  OF  A  GAS-VAPOR  BUBBLE 

IN  A  LIQUID-GAS  SOLUTION  SUBJECT  TO  A  CENTRI¬ 
FUGAL  FIELD .  17 

2.1  The  Equilibrium  Conditions  .  17 

2.2  The  Equilibrium  Shape  of  the  Gas-Vapor 

Bubble .  35 

3  THE  NATURE  OF  THE  EQUILIBRIUM  STATE .  4  7 

4  THE  CRITICAL  SIZE  OF  A  GAS-VAPOR  BUBBLE  IN  A 

DILUTE  SOLUTION  AT  NEGATIVE  PRESSURES  ...  73 

5  THE  CRITICAL  RADIUS  OF  A  GAS-VAPOR  BUBBLE  IN  A 

DILUTE  SOLUTION  OF  AIR  IN  WATER  AT  NEGATIVE 
PRESSURES .  8  9 

6  SUMMARY  AND  CONCLUSIONS  .  116 

BIBLIOGRAPHY  .  122 

APPENDIX  A:  THE  GENERAL  EXPRESSION  FOR  THE  EQUILIBRIUM 

SHAPE  OF  THE  GAS-VAPOR  BUBBLE .  12  9 

B:  DERIVATION  OF  EXPRESSIONS  FOR  THE  SURFACE 

AREA,  VOLUME,  AND  MOMENT  OF  INERTIA  OF  THE 

SPHEROID .  133 

C:  NUMERICAL  DATA  USED  IN  CALCULATIONS  .  .  .  140 


Vll 


' 


LIST  OF  TABLES 


TABLE  DESCRIPTION 

1  The  Dependence  of  the  Critical  Radius  on 
the  Air  Content  for  Water  with  Various 
Negative  Pressures  at  a  Temperature  of  25°C 

2  The  Temperature  Variation  of  the  Surface 

Tension,  Vapor  Pressure,  and  Henry's  Law 
Constant  for  a  Dilute  Solution  of  Air  in 
Water  . 

3a  The  Percentage  Reductions  in  the  Partial 
Pressures  for  Various  Negative  Pressures 

3b  The  Percentage  Reductions  in  the  Partial 
Pressures  for  Various  Concentrations  .  . 

3c  The  Percentage  Reductions  in  the  Partial 
Pressures  for  Various  Temperatures  .  .  . 


viii 


PAGE 

105 

106 

112 

113 

114 


LIST  OF  FIGURES 


FIGURE  PAGE 

1-1  Density  Fluctuations  in  a  Liquid  with  Time  4 

1-2  Galloway's  Experimental  Results  for  the 

Dependence  of  the  Cavitation  Threshold 
Pressure  on  the  Percentage  Gas  Concentration  10 

1-3  Results  of  Donoghue  et:  al^  for  the  Dependence 

of  the  Rupture  Strength  of  Benzene  on  the 
Dissolved  Gas  Content  .  12 

2.1- 1  The  Thermodynamic  System .  18 

2.1- 2  The  Three  Major  Subsystems  .  21 

2.2- 1  A  Plane  Section  Through  the  Axis  of  Rotation 

of  the  Gas-Vapor  Bubble .  3  6 

3-1  A  Comparison  of  the  Plane  Sections  of  the 

Gas-Vapor  Bubble  in  the  Equilibrium  and  Non- 
Equilibrium  States .  48 

3- 2  The  Behavior  of  the  Helmholtz  Potential  in  the 

Neighborhood  of  the  Equilibrium  State  for 
Which  e  Equals  Zero .  66 

4- 1  A  Comparison  of  the  Primary  and  Reference 

Systems .  8  0 

5- 1  Deviation  of  the  Gas-Vapor  Bubble  From 

Sphericity .  91 

5-2a  Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  the  Reference  Pressure  at 
Various  Negative  Pressures.  For  a  solution 
of  air  in  water  at  a  concentration  of  10-4and 
a  temperature  of  25°C .  93 

5-2b  Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  the  Reference  Pressure  at 
Various  Concentrations.  For  a  solution  of  air 
in  water  at  a  negative  pressure,  P^,of  -10.13 
MPa  and  a  temperature  of  25°C .  94 


IX 


FIGURE 


PAGE 


5-2c 


5-3a 


5-3b 


5-3c 


5-3d 


5-4 


5-5a 


5-5b 


Theoretical  Results  for  the  Dependence  of 
the  Critical  Radius  on  the  Reference  Pressure 
at  Various  Temperatures.  For  a  solution  of 
air  in  water  at  a  negative  pressure,  Pr; ,  of 
-10.13  MPa  and  a  concentration  of  10-^  ...  95 

Theoretical  Results  for  the  Dependence  of 
the  Critical  Radius  on  Negative  Pressure 
at  Various  Concentrations.  For  a  solution  of 
air  in  water  at  a  temperature  of  25°C  ...  98 

Theoretical  Results  for  the  Dependence  of  the 

Critical  Radius  on  Negative  Pressure  at 

Various  Concentrations.  For  a  solution  of 

air  in  water  at  a  temperature  of  100°C  ...  99 

Theoretical  Results  for  the  Dependence  of  the 

Critical  Radius  on  Negative  Pressure  at 

Various  Concentrations.  For  a  solution  of 

air  in  water  at  a  temperature  of  200°C  .  .  .  100 

Theoretical  Results  for  the  Dependence  of  the 

Critical  Radius  on  Negative  Pressure  at 

Various  Concentrations.  For  a  solution  of 

air  in  water  at  a  temperature  of  300°C  .  .  .  101 

Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Dissolved  Gas  Content  at 
Various  Negative  Pressures.  For  a  solution 
of  air  in  water  at  a  temperature  of  25°C  .  .  103 

Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Temperature  for  a  Pure 
Vapor  Bubble  in  Water  at  Various  Negative 
Pressures .  107 

Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Temperature  at  Various 
Negative  Pressures.  For  a  solution  of  air  in 
water  at  a  concentration  of  10”^ .  108 


x 


. 


LIST  OF  SYMBOLS 


A  surface  area 

a  maximum  radius  of  the  gas-vapor  bubble 

B  eccentricity  of  an  ellipse 

b  half-length  of  the  gas-vapor  bubble 

C  concentration  of  gas  in  a  liquid 

Cq  reference  concentration  of  gas  in  a  liquid 

Cg  equilibrium  concentration  of  gas  in  the  liquid 

in  the  reference  system 

C  concentration  of  the  solution  at  the  membrane  in 

SO  the  reference  system 

D  elemental  quantity 

D  diffusion  coefficient 

c 

e  shape  parameter 

F  Helmholtz  potential 

f  dimensionless  parameter 

G  Gibb's  potential 

H2  radius  of  curvature  corresponding  to  the  principal 

curvature 

h  radius  of  curvature  at  the  center  of  the  gas-vapor 

bubble 

I  moment  of  inertia  of  a  body  about  a  cylindrical  axis 

K  modulus  of  elliptic  integral 

K  Henry's  law  constant 

H 

k  Boltzmann's  constant 

l  half-length  of  the  centrifugal  system 


xi 


molecular  weight  of  component  i 
mass 

moles  of  component  i 

total  number  of  components 

normal  component  of  motion  of  the  surface 
element  DA 

pressure 

pressure  in  the  liquid-gas  solution  at  the  center 
of  the  system 

partial  pressure  of  component  i  in  the  gas-vapor 
mixture 

reference  pressure  for  the  centrifugal  system 
probability 

general  reference  pressure 

vapor  pressure  of  the  solvent 

dimensionless  radial  co-ordinate 

critical  radius  of  a  spherical  gas-vapor  bubble 

radial  co-ordinate 

entropy 

arc  length  of  a  curve 
temperature 
internal  energy 
volume 

partial  molar  volume  of  the  solute  in  a  solution 
specific  volume  of  the  solvent 
reversible  work 


dimensionless  axial  co-ordinate 


mole  fraction  of  component  i 
axial  co-ordinate 
transformation  variable 
transformation  variable 


adsorption  of  component  i 

surface  tension  of  the  gas-vapor  bubble 

virtual  variation 

difference  in  the  maximum  radius  between  the 
equilibrium  and  non-equilibrium  systems 

dimensionless  variables 
angular  co-ordinate 

angle  associated  with  the  slope  of  a  curve 
principal  curvatures  of  the  "dividing  surface" 
undetermined  multiplier 
chemical  potential  of  component  i 

dimensionless  variables 

density 

characteristic  time  scale  for  diffusion 

angle  between  the  normal  to  the  curve  and  the 
positive  x-axis 

angle  between  the  normal  and  radial  directions 
angle  associated  with  elliptic  integral 
undetermined  function  of  P  and  T 
angular  velocity  of  the  system 


Subscripts 


i  component  i 

j  component  j 

a  equilibrium  state 

3  non-equilibrium  state 

v  equilibrium  state  without  the  gas-vapor 

bubble  present 

Superscripts 

G  gas-vapor  mixture 

GR  gas-vapor  mixture  in  the  reference  system 

L  liquid-gas  solution 

LR  liquid-gas  solution  in  the  reference  system 

S  "dividing  surface" 

T  total  quantity 


xiv 


CHAPTER  1 


INTRODUCTION 

The  transition  of  a  fluid  from  a  liquid  phase  to  a 
vapor  phase  characterizes  many  important  fluid  phenomena. 

Two  of  the  most  interesting  of  these  from  an  engineering 
point  of  view  are  boiling  and  cavitation.  Boiling  and 
cavitation  can  be  described  as  the  explosive  vaporization 
of  the  liquid  into  a  vapor-filled  cavity  [1] .  Boiling  is 
caused  by  heating  under  constant  pressure,  while  cavitation 
is  the  result  of  a  reduction  in  pressure  at  essentially 
constant  temperature.  The  onset  of  both  of  these  phase 
transition  phenomena  is  not  as  yet  fully  understood.  For 
example,  cavitation  is  said  to  occur  when  the  pressure  in 
the  fluid  drops  below  the  vapor  pressure  of  the  liquid.  In 
actual  practice,  this  criterion  is  much  too  simple  to  be 
able  to  accurately  predict  the  onset  of  cavitation. 

Boiling  and  cavitation  represent  phase  transitions 
as  a  result  of  a  metastable  state  in  a  liquid.  A  metastable 
state  is  stable  for  infinitesimal  deviations  about  the 
equilibrium  state;  however,  for  a  finite  deviation,  the  system 
does  not  return  to  the  original  state,  but  proceeds  to  some 
other  equilibrium  state  for  which  the  entropy  of  the  system 
is  greater  than  in  the  original  state.  Given  enough  time,  a 
system  will  move  from  a  metastable  state  to  a  stable  one. 
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The  mechanism  of  this  phase  transition  is  crucial  to  an 
understanding  of  boiling  and  cavitation. 

Both  boiling  and  cavitation  require  the  formation 
of  a  vapor  cavity  of  particular  size  in  the  metastable  liquid 
state.  According  to  Plesset  [2] ,"a  central  problem  in 
cavitation  and  boiling  is  how  macroscopic  vapor  cavities 
can  form  when  moderate  tensions  are  applied  to  a  liquid".  If 
the  liquid  pressure  is  mechanically  reduced  sufficiently 
below  the  vapor  pressure,  the  liquid  can  enter  a  state  of 
tension,  or  negative  pressure.  Negative  pressures  are  made 
possible  in  a  liquid  by  intermolecular  attraction.  In  this 
sense,  a  liquid  more  closely  resembles  an  amorphous  solid 
than  a  gas. 

Theoretical  calculations  have  shown  that  the  fracture 
of  a  pure  liquid  leading  to  the  creation  of  a  vapor  cavity 
could  only  be  produced  by  extremely  large  tensions^"  or  super¬ 
heats  [4].  However,  experimentally  much  lower  tensile 
strengths  and  superheat  limits  have  been  obtained  for  common 
engineering  fluids.  In  an  attempt  to  explain  the  relatively 
low  tensions  and  superheats  commonly  encountered,  the  presence 
of  nuclei  has  been  postulated  in  ordinary  liquids.  These 
nuclei  would  serve  to  initiate  the  formation  of  the  vapor 
cavity . 

It  is  well  known  that  a  "dirty"  liquid  has  a  much 
lower  tensile  strength  than  the  pure  liquid.  Heterogeneous 

■'"For  water  at  300°K,  Fischer  [3]  calculated  a  tensile 
limit  of  approximately  134  MPa. 
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nucleation  theory  considers  a  foreign  particle  or  surface 
to  serve  as  the  nucleation  site.  Thus,  the  presence  of  a 
surface  or  colloidal  impurity  in  the  liquid  would  be 
responsible  for  the  reduction  in  tensile  strength. 

Another  possible  reason  for  the  low  liquid  tensions 
experimentally  obtained  is  the  presence  of  dissolved  gas. 
Even  when  careful  attention  is  given  to  the  purity  of  the 
liquid,  in  practice  small  amounts  of  dissolved  gas  are 
usually  present.  To  date,  the  effect  of  the  presence  of 
dissolved  gas  on  the  tensile  strength  of  a  liquid  has  not 
been  conclusively  determined.  There  is  general  disagreement 
as  to  whether  the  gas  significantly  lowers  the  tensile 
limit  [5] . 

The  presence  of  dissolved  gas  is  important  when  the 
vapor  cavity  formation  is  initiated  by  homogeneous  bubble 
nucleation.  Homogeneous  nucleation  theory  [6,7]  postulates 
that  a  thermodynamic  fluctuation  in  a  metastable  state  pro¬ 
duces  an  "embryo"  of  the  new  phase,  which  if  sufficiently 
large  in  size  serves  to  initiate  the  phase  transition.  For 
the  transition  of  a  liquid  to  a  vapor,  the  embryo  can  be 
regarded  as  a  very  small  vapor  bubble  or  nucleus.  From 
statistical  physics  [8] ,  it  is  understood  that  the  instan¬ 
taneous  value  of  the  density  of  a  liquid  fluctuates  over  a 
range  of  values.  Figure  1-1  illustrates  the  density  at  a 
particular  location  in  a  liquid  as  a  function  of  time. 
Although  the  instantaneous  value  fluctuates  rapidly,  the 
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Figure  1-1  Density  Fluctuations  in  a  Liquid  with  Time 
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time  average  value  remains  constant  for  a  stable  equilibrium 
state.  Small  deviations  in  the  density  do  not  disturb  the 
equilibrium  state  of  the  liquid.  However,  a  sufficiently 
large  fluctuation  in  the  density  of  the  liquid  can  cause  a 
vapor  nucleus  to  be  formed.  For  a  liquid  phase  in  a  meta¬ 
stable  state,  it  is  possible  that  the  nucleus  will  not  decay 
as  it  does  for  a  stable  equilibrium  state,  but  instead  grow 
to  form  a  center  for  vaporization. 

Whether  the  nucleus  grows  or  decays  is  explained  by 
the  concept  of  a  critical  size,  or  in  the  case  of  a  spherical 
nucleus,  a  critical  radius.  There  is  a  critical  size  for 
which  a  nucleus  in  a  metastable  state  will  grow  to  become  a 
center  for  the  formation  of  a  new  phase.  If  the  nucleus  is 
less  than  the  critical  size,  the  nucleus  will  decay.  How¬ 
ever,  if  the  nucleus  is  greater  than  the  critical  size,  the 
nucleus  will  grow  to  become  in  this  case  a  center  for  vapor¬ 
ization  . 

From  statistical  physics  [9],  it  is  possible  to 
predict  the  probability  that  a  critical  size  nucleus  will 
occur.  The  probability,  Pr,  is  inversely  proportional  to  the 
exponential  of  the  reversible  work,  W  ,  required  to  form  the 
nucleus,  i.e. 

WR 

Pr  “  exp  (-  (1-1) 

where  k  is  Boltzmann's  constant,  and  T  is  the  temperature. 

The  reversible  work  in  the  above  expression  is  a  function  of 
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the  critical  size  of  the  vapor  nucleus. 

Although  the  discussion  has  referred  so  far  to  a 
pure  liquid,  it  is  of  special  interest  to  consider  the  case 
of  a  liquid-gas  solution  in  order  to  include  the  effect  of 
the  dissolved  gas  on  the  system.  The  reversible  work 
required  to  form  a  critical  size  nucleus  in  a  liquid-gas 
solution  can  be  determined  if  the  size  of  the  nucleus  has 
been  obtained  as  a  function  of  the  properties  of  the  solution. 
For  a  spherical  nucleus  of  radius  R  ,  in  a  two  component 
solution,  a  useful  representation  is 

R  =  R  (T,P,C)  (1-2) 

where  T,  P,  and  C  are,  respectively,  the  temperature, 
pressure,  and  concentration  of  the  solution. 

Ward,  Balakrishnan ,  and  Hooper  [10]  have  developed 
an  expression  for  the  critical  radius  of  a  gas-vapor  bubble 
in  a  liquid-gas  solution  at  positive  pressure.  The  condi¬ 
tions  of  thermodynamic  equilibrium  for  the  general  case  of  a 
multicomponent  solution  yielded  the  equality  of  chemical 
potentials,  and  the  Laplace  equation  of  capillarity  relating 
the  pressure  difference  across  the  interface  to  the  curvature 
of  the  bubble  at  that  point.  For  the  special  case  of  a  dilute 
solution  of  a  gas  in  a  liquid,  the  critical  radius  of  a 
spherical  bubble  was  expressed  as  a  function  of  the  tempera¬ 
ture,  concentration,  and  pressure  of  the  dilute 
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solution,  i.e. 


where 


R 

c 


P 


n  =  exp  (p_pv)  " 


(1-2) 


(1-3) 


The  surface  tension  y ,  specific  volume  vL,  and  vapor 
pressure  are  functions  of  the  temperature,  T.  The  refer¬ 
ence  concentration  Cq  is  that  of  the  liquid  saturated  with 
the  gas  across  a  flat,  non-rigid  membrane  permeable  only  to 
the  gas  component.  The  reference  system  is  strictly  limited 
to  positive  pressures  for  the  liquid.  If  the  liquid  is  in  a 
state  of  tension,  the  equality  of  pressure  across  the  membrane 
would  require  the  gas  to  also  exist  at  a  negative  pressure, 
which  is  impossible.  Thus,  the  expression  for  the  critical 
radius  given  by  equation  (1-2)  is  only  valid  for  positive 
pressures  in  the  solution. 

If  homogeneous  nucleation  theory  is  to  be  used  to 
explain  the  onset  of  cavitation  and  boiling  in  liquid-gas 
solutions,  the  expression  for  the  critical  size  of  a  gas-vapor 
nucleus  must  be  extended  to  include  the  possibility  of 
negative  pressures.  From  this  expression,  the  reversible  work 
required  to  form  a  critical  size  nucleus  in  such  a  metastable 

2 

Since  Ward  et  alL.  included  the  effect  of  non-ideal 
gas  behavior  inside  the  bubble,  the  expression  obtained  for 
R  included  the  activity  coefficients  of  each  component.  Ex¬ 
pression  (1-2)  results  from  setting  the  activity  coefficients 
equal  to  one  in  the  original  expression,  i.e.  it  is  assumed  that 
the  gas  and  vapor  inside  the  bubble  behave  as  an  ideal  gas 
mixture . 
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state  can  be  calculated.  Once  the  reversible  work  is  known, 
the  probability  of  the  formation  of  a  critical  size  nucleus 
can  be  obtained.  This  investigation  will  consider  an  ex¬ 
pression  for  the  critical  size  of  a  gas-vapor  nucleus  in  a 
liquid-gas  solution  at  negative  pressure,  and  relate  it  to 
the  reversible  work  required  to  form  the  nucleus. 

Since  this  investigation  will  involve  a  solution  at 
negative  pressure,  it  is  useful  to  consider  some  of  the 
techniques  that  have  been  developed  for  producing  tensions 
in  liquids.  The  principal  purpose  of  experimentally  gener¬ 
ating  negative  pressures  in  liquids  has  been  the  measurement 
of  liquid  tensile  strengths.  A  thorough  review  of  the  methods 
used  to  determine  the  tensile  strengths  of  liquids,  and  some 
of  the  results  obtained  are  given  by  Blake  [11] . 

One  method  of  generating  liquid  tensions  is  by  the 
use  of  pressure  fields.  For  example,  the  region  of  rare¬ 
faction  for  an  acoustical  standing  wave  in  a  liquid  can 
exhibit  negative  pressures.  Galloway  [12]  used  this  technique 
to  study  sonically  induced  cavitation  in  water  and  benzene, 
among  other  liquids.  Radially  symmetric  standing  waves 
produced  by  a  spherical  resonator  were  used  to  generate 
acoustic  pressures  in  the  body  of  the  liquid  away  from  any 
solid  surface.  Dynamic  tensile  strengths  of  approximately 
20  MPa  for  water  and  14  MPa  for  benzene  were  obtained.  It 
was  determined  that  the  air  content  of  the  liquid  signifi¬ 
cantly  limited  the  maximum  negative  pressure  the  liquid  could 
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withstand.  Figure  1-2  illustrates  Galloway's  results  for 
the  dependence  of  the  cavitation  threshold  pressure  of  water 
and  benzene  on  the  percentage  air  concentration. 

A  centrifugal  pressure  field  can  also  be  used  to 
generate  liquid  tensions.  A  description  of  the  method  used 
by  Briggs  [13]  will  suffice  to  explain  the  technique  involved. 
The  liquid  was  contained  in  a  thin  capillary  tube  open  at 
both  ends.  About  one  centimeter  from  each  end,  the  tube  was 
bent  back  upon  itself  through  an  angle  of  140°  to  form  a  Z 
shape.  The  tube  was  mounted  on  a  spinner  connected  to  the 
end  of  a  vertical  shaft  of  an  electric  motor.  The  entire 
system  was  enclosed  in  a  large  cylinder.  When  operating, 
the  air  pressure  inside  the  cylinder  was  reduced  to  4  or  5 
kPa  (3  or  4  cm  Hg) .  The  tube  was  mounted  on  a  horizontal 
plane  with  the  center  intersecting  the  projected  spin  axis. 
Thus,  when  spinning,  each  half  column  of  liquid  pulled 
against  the  other,  placing  the  liquid  at  the  center  in  a 
state  of  tension.  The  fracture  of  the  column  of  liquid  at 
the  mid-point  of  the  tube  indicated  that  the  maximum  tensile 
strength  had  been  reached.  Knowledge  of  the  length  of  the 
tube  and  the  angular  velocity  of  the  system  enabled  the 
maximum  negative  pressure  to  be  calculated  for  a  liquid  of 
given  density.  Briggs  notes  that  "scrupulous  cleanliness" 
was  necessary. 

Reynolds  [14] ,  in  1878,  was  the  first  to  use  a  cen¬ 
trifugal  system  to  measure  the  tensile  strength  of  a  liquid. 


Cavitation  Threshold  Pressure  (MPa) 
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Galloway's  Experimental  Results  for  the 
Dependence  of  the  Cavitation  Threshold  Pressure 
on  the  Percentage  Gas  Concentration.  (100% 
corresponds  to  the  saturation  concentration  at 
a  pressure  of  0.10  MPa  and  a  temperature  of  22°C) 


Figure  1-2 
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He  obtained  a  maximum  tension  of  about  0.49  MPa  for  ordinary 
tap  water.  Temperley  and  Chambers  [15] ,  in  an  attempt  to 
repeat  his  results  in  1946,  obtained  a  maximum  value  of 
0.57  MPa  for  water  at  room  temperature,  The  largest  experi¬ 
mentally  measured  tension  for  water  is  that  due  to  Briggs. 

Using  the  centrifugal  system  described  above,  he  obtained  a 
maximum  tension  of  about  27.7  MPa  for  water  around  10 °C. 

A  centrifugal  system  has  also  been  applied  to  the 
study  of  liquids  other  than  water.  Apfel  and  Smith  [16] 
obtained  a  maximum  tension  of  17.6  MPa  for  di-ethyl  ether 
using  an  air-driven  centrifugal  device.  The  significance 
of  this  value  is  that  it  lies  "between  88  and  93  per  cent  of 
the  tensile  strength  of  di-ethyl  ether  at  room  temperature 
as  predicted  by  homogeneous  nucleation  theory."  Donoghue , 
Vollrath,  and  Gerjuoy  [17]  obtained  a  maximum  tension  of 
15.9  MPa  for  benzene  using  a  system  similar  to  that  of  Briggs. 
Like  Galloway,  Donoghue  ejt  al.  observed  the  effect  of  dissolved 
gas  on  the  tensile  strength  of  benzene.  Some  of  their 
results  are  presented  in  Figure  1-3.  They  indicate  that  the 
rupture  strength  is  decreased  by  the  addition  of  dissolved 
gas . 

Although  the  results  of  Galloway  and  Donoghue  et  al . 
suggest  that  the  presence  of  dissolved  gas  decreases  the 
tensile  strength  of  a  liquid,  their  results  are  not  conclu¬ 
sive.  In  a  brief  discussion  of  Galloway's  results,  Hammitt 
[18]  notes  that  for  a  higher  gas  content,  i.e.  fifty  to 
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Tension  at  Rupture  (MPa) 


12 


Log  Gas  Pressure  (mm  Hg) 


Results  of  Donoghue  et  al.  for  the  Dependence 
of  the  Rupture  Strength  of  Benzene  on  the 
Dissolved  Gas  Content 


Figure  1-3 
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ninety-five  per  cent  saturation,  the  results  of  Iyengar  and 
Richardson  [19]  using  ultrasonic  cavitation  in  static  systems 
show  no  effect  of  dissolved  air  on  the  cavitation  threshold. 
Kuper  and  Trevena  [20]  reported  on  theoretical  grounds  that 
"the  presence  of  dissolved  gas  only  causes  a  negligible 
diminution  in  the  intrinsic  strength  of  liquids."  Harvey 
[21]  found  that  prepressurizing  the  liquid  to  drive  the  free 
gas  into  solution  enabled  considerable  tensions  to  be  obtain¬ 
ed  in  water.  He  demonstrated  that  in  considering  the  effect 
of  gas  content  on  the  tensile  strength  of  liquids,  one  must 
discriminate  between  dissolved  and  undissolved  gas. 

In  general,  it  should  be  noted  that  the  use  of  a 
centrifugal  system  to  measure  the  tensile  strength  of  a 
liquid  requires  that  the  rupture  originate  in  the  body  of  the 
liquid.  However,  in  many  cases  it  has  not  been  clearly 
established  whether  the  liquid  column  fails  due  to  loss  of 
adhesion  at  the  walls  of  the  tube  or  loss  of  cohesion  in  the 
body  of  the  liquid. 

This  investigation  will  consider  the  use  of  a  centri¬ 
fugal  field  to  produce  a  region  of  negative  pressure  in  a 
liquid.  If  a  gas-vapor  nucleus  is  imagined  to  exist  at  the 
center  of  rotation  of  a  system  similar  to  that  of  Briggs, 
the  solution  outside  the  nucleus  would  be  subject  to  negative 
pressures.  Thus,  the  thermodynamic  analysis  of  a  gas-vapor 
nucleus  located  at  the  axis  of  rotation  of  a  narrow  tube 
filled  with  a  liquid-gas  solution  will  yield  an  expression 
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for  the  critical  size  similar  to  that  obtained  by  Ward  et  al ., 
valid  for  the  case  of  negative  pressure.  The  expression 
for  the  critical  size  of  the  nucleus  in  terms  of  the  proper¬ 
ties  of  the  solution  will  require  the  use  of  a  special 
reference  system  different  from  that  of  Ward  et  ad. ,  which 
is  valid  for  a  solution  at  negative  pressure. 

For  a  system  subject  to  a  centrifugal  field,  the 
deformation  of  the  bubble  must  be  considered.  A  gas-vapor 
bubble  rotating  about  a  central  axis  cannot  be  assumed  to 
have  a  spherical  shape.  Rather,  the  effect  of  the  centrifugal 
field  would  be  such  as  to  produce  a  deformation  consistent 
with  the  axial  symmetry  of  the  system.  Lord  Rayleigh  [22] 
considered  the  effect  of  "capillary  tension"  on  an  incom¬ 
pressible  fluid  mass  rotating  with  uniform  angular  velocity. 

He  solved  the  differential  equation  for  the  system  using  a 
power  series  expansion.  When  only  first  order  terms  were 
considered,  the  shape  of  a  plane  section  through  the  axis  of 
rotation  of  the  bubble  was  approximated  by  an  ellipse. 

Rosenthal  [23]  and  Princen,  Zia,  and  Mason  [24]  both 
considered  the  shape  of  a  fluid  bubble  immersed  in  and 
rotating  with  a  heavier  fluid.  Beginning  with  Laplace's 
equation  of  capillarity,  both  obtained  solutions  for  the  shape 
of  a  plane  section  through  the  axis  of  rotation  in  terms  of 
elliptic  integrals  of  the  first  and  second  kind.  Each  method 
differed  in  the  choice  of  the  bubble  parameter  to  which  the 
size  and  shape  were  referenced,  although  the  results  were 


, 
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easily  related.  Both  concluded  that  for  high  angular  veloc¬ 
ities,  the  shape  of  the  bubble  became  essentially  that  of  a 
cylinder  with  rounded  ends. 

Princen  et  ad.  used  their  expression  for  the  bubble 
length  to  determine  the  interfacial  tension  of  an  experiment¬ 
al  system.  Measuring  the  length  of  a  bubble  of  known  volume 
at  a  given  speed  enabled  them  to  calculate  the  surface  tension. 
In  the  course  of  their  experiments,  they  observed  the  bubble 
to  become  cylindrical  in  shape  for  large  angular  velocities, 
as  they  had  predicted. 

The  conclusions  of  Rosenthal  and  Princen  et  al.  are 
not  directly  applicable  since  their  systems  treated  fluid 
bubbles  of  constant  volume.  That  is,  a  bubble  of  given 
volume  initially  placed  in  their  systems  at  rest,  will  retain 
that  volume  even  when  accelerated  to  different  angular 
velocities.  For  a  gas-vapor  bubble  in  a  liquid-gas  solution, 
the  effect  of  mass  transfer  across  the  interface  must  be 
included.  The  equilibrium  conditions  will  determine  the 
distribution  of  mass  inside  and  outside  the  bubble. 

In  this  investigation,  the  conditions  of  thermo¬ 
dynamic  equilibrium  will  be  obtained  for  a  gas-vapor  nucleus 
at  the  center  of  rotation  of  a  centrifugal  system.  A  stability 
analysis  will  determine  the  nature  of  the  equilibrium  state. 
From  the  conditions  of  equilibrium,  the  shape  of  the  bubble 
will  be  obtained  for  negative  pressures  in  the  solution  out¬ 
side  the  bubble.  For  the  specific  case  of  a  gas-vapor  bubble 
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in  a  dilute  solution,  the  use  of  an  appropriate  reference 
condition  will  enable  the  critical  size  to  be  expressed  in 
terms  of  the  properties  of  the  dilute  solution. 


CHAPTER  2 


THE  EQUILIBRIUM  STATE  OF  A  GAS-VAPOR  BUBBLE 
IN  A  LIQUID-GAS  SOLUTION  SUBJECT  TO 
A  CENTRIFUGAL  FIELD 

2 . 1  The  Equilibrium  Conditions 

This  section  considers  the  conditions  of  thermo¬ 
dynamic  equilibrium  for  a  gas-vapor  bubble  in  a  multi- 
component  liquid-gas  solution  at  negative  pressure.  In 
order  to  place  the  liquid  phase  outside  the  bubble  in  a 
state  of  tension,  the  entire  system  is  considered  subject 
to  a  centrifugal  field. 

The  thermodynamic  system  is  represented  by  the 
schematic  in  Figure  2.1-1.  A  long,  narrow  tube  rotates  at 
constant  velocity  around  a  central  axis  denoted  x.  The  tube 
is  filled  with  a  solution  of  n-1  gases  in  an  incompressible 
liquid.  Each  chemical  component  is  considered  to  be  inert, 
i.e.  there  are  no  chemical  reactions  in  the  system.  The 
pressure  in  the  liquid  solution  at  both  ends  of  the  tube — a 
radial  distance  £from  the  x-axis — is  denoted  by  the  reference 
pressure  P  .  A  single  gas-vapor  bubble  is  located  at  the 
center  of  the  tube  in  the  liquid-gas  solution.  The  liquid- 
gas  solution  containing  the  gas-vapor  bubble  composes  the 
thermodynamic  system  considered.  The  entire  system  is  assumed 
to  rotate  with  angular  velocity  w,  and  there  are  no  other 
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Figure  2.1-1  The  Thermodynamic  System 
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fields  acting  on  the  system.  For  the  closed,  constant- 
volume  system  considered,  the  total  mass  and  total  volume 
remain  constant.  The  entire  system  is  in  diathermal  contact 
with  a  heat  reservoir  at  temperature  T. 

The  equilibrium  state  of  a  thermodynamic  system  is 

that  which  minimizes  the  value  of  the  internal  energy  for  a 

fixed  value  of  the  total  entropy.  However,  for  a  system 

at  constant  temperature,  it  is  more  convenient  to  use  the 

Helmholtz  potential.^"  The  natural  variables  for  the 

Helmholtz  potential,  F,  are  the  temperature,  T,  volume,  V, 

2 

and  mole  numbers,  N.  The  equilibrium  state  for  a  system  at 
constant  temperature  is  that  which  minimizes  the  value  of  the 
Helmholtz  potential  at  that  temperature  [25].  The  advantage 
of  using  the  Helmholtz  potential  formulation  is  a  reduction 
by  one  in  the  number  of  independent  variables  due  to  the 
constant  temperature  condition.  Given  an  expression  for  the 
total  free  energy  of  the  system,  the  equilibrium  conditions 
are  obtained  by  setting  the  variation  of  the  energy  at 
constant  temperature  equal  to  zero. 

^The  Helmholtz  potential  is  the  Legendre  trans¬ 
formation  of  the  internal  energy,  U,  replacing  the  entropy, 

S,  by  the  temperature,  T,  as  an  independent  variable,  i.e. 

F  =  U  -  TS 

2 

The  number  of  molecules,  expressed  in  moles,  of  a 
chemical  component  i  is  denoted  by  the  mole  number  N..  The 
notation  N  is  used  to  represent  the  mole  numbers  of  each  of 
the  respective  components,  i.e. 


Nf,  N2, 


.  .  .  , 


N  . 

n 
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The  first  step  in  the  equilibrium  analysis  is  to 
formulate  an  expression  for  the  total  free  energy  of  the 
system  using  the  Helmholtz  potential  representation.  The 
total  system  represented  in  Figure  2.1-1  is  comprised  of 
three  subsystems:  the  liquid-gas  solution,  the  gas-vapor 
mixture,  and  the  interface,  denoted  by  the  superscripts  L, 

G,  and  S,  respectively.  Figure  2.1-2  illustrates  the  three 
subsystems  created  by  the  bubble  surface.  The  properties  of 
a  single  subsystem  are  not  uniform  and  homogeneous  through¬ 
out  due  to  the  influence  of  the  centrifugal  field.  In  order 
to  account  for  the  spatial  variation  of  the  local  properties, 
small  elements  of  each  subsystem  are  considered.  The  symbol 

D  is  used  to  indicate  small  quantities  of  the  extensive 
3 

variables  for  which  the  properties  are  constant. 

The  Helmholtz  potential  in  the  liquid-gas  solution 
and  the  gas-vapor  mixture  can  be  expressed  as  a  function  of 
the  local  values  of  the  temperature,  volume,  and  mole  numbers 
of  the  chemical  components.  The  expression  for  the  Helmholtz 
potential  of  a  small  element  of  either  subsystem  at  constant 
temperature  is 

n 

DF  =  -  PDV  +  Z  y  .  DN .  (2.1-1) 

i=l  1  1 

where  P  is  the  pressure,  and  y^  is  the  chemical  potential  of 
the  ith  component.  Equation  (2.1-1)  can  be  integrated  over 
each  subsystem  to  obtain  an  expression  for  the  Helmholtz 

An  extensive  variable  or  parameter  is  one  that  is 
additive  over  the  extent  of  the  system,  i.e.  V  or  ISh. 
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Figure  2.1-2  The  Three  Major  Subsystems 
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potential  of  the  liquid-gas  solution  and  the  gas-vapor 
mixture,  i.e. 


and 


PLDVL  + 


n 

E 

i=l 


PGDVG  +  E 
i=l 


(2.1-2) 


(2.1-3) 


The  expression  for  the  Helmholtz  potential  associated 

with  the  interface  is  obtained  from  surface  thermodynamics 

[26] .  The  interface  is  modelled  as  a  mathematical  surface  of 

discontinuity  called  the  "dividing  surface".  The  position  of 

the  "dividing  surface"  is  so  chosen  as  to  eliminate  the 

dependence  of  the  free  energy  associated  with  the  interface 

on  the  curvature  of  the  surface  [27].  It  has  zero  thickness, 

and  hence  no  volume.  The  area  replaces  the  volume  as  an 

independent  parameter  in  the  expression  for  the  Helmholtz 

4 

potential.  Excess  quantities  are  associated  with  the 
"dividing  surface"  in  an  analogous  way  to  their  counterparts 
in  the  other  two  systems.  Thus,  the  Helmholtz  potential  in 
the  interface  is  expressed  as  a  function  of  the  temperature, 
area,  and  excess  mole  numbers  of  the  chemical  components. 

^An  excess  quantity  represents  the  difference  between 
the  actual  amount  of  an  extensive  quantity  contained  in  the 
interfacial  region,  and  the  amount  that  would  be  contained  if 
the  phases  on  either  side  of  the  "dividing  surface"  were 
assumed  to  extend  right  up  to  the  surface  in  a  continuous 
and  uniform  manner. 
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The  expression  for  the  Helmholtz  potential  of  a  small  element/ 
DA,  of  the  "dividing  surface"  at  constant  temperature  is 

DFS  =  yDA  +  E  y.S  DN.S  (2.1-4) 

.-,11 

i=l 


where  y  is  the  surface  tension  [28] .  Integrating  over  the 
"dividing  surface"  gives  the  expression  for  the  Helmholtz 
potential  of  the  interface,  i.e. 


yDA  + 
A 


n 

I 

i=l 


(2.1-5) 


The  Helmholtz  potential  of  the  entire  system  is 
simply  the  sum  of  the  potentials  of  the  three  subsystems, 

i.e. 

F  =  FL  +  FG  +  FS  (2.1-6) 

Using  equations  (2.1-2)  ,  (2.1-3)  ,  and  (2.1-5)  in  (2.1-6) 

gives  an  expression  for  the  intrinsic  free  energy  of  the 
total  system,  i.e. 
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F  = 


pldvl  + 
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i=l 
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L  L 

y .  dn. 

1  1 


N.' 
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(2.1-7) 


The  system  shown  in  Figure  2.1-1  is  subject  to  a 
centrifugal  field  generated  by  a  rotation  about  the  central 
axis.  The  expression  for  the  intrinsic  free  energy  given 
by  equation  (2.1-7)  does  not  include  the  energy  contribution 
due  to  the  field.  The  potential  energy  associated  with  a 
small  element  of  mass.  Dm,  in  the  presence  of  the  field  is 

-  -j  oo2r2  Dm  (2.1-8) 

where  u)  is  the  angular  velocity,  and  r  is  the  radial 

5 

distance  of  element  Dm  from  the  axis  of  rotation. 

- 

The  method  whereby  the  effect  of  the  centrifugal 
field  is  introduced  closely  follows  that  developed  by 
Gibbs  [29]  in  treating  the  equilibrium  of  a  system  under 
the  influence  of  gravity. 
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Integrating  expression  (2.1-8)  over  each  subsystem  gives 
the  contribution  of  the  centrifugal  field  to  the  free 
energy  of  each  subsystem,  i.e. 

-  j  co2r2  Dm  (2.1-9) 

.  m 

Finally,  combining  expressions  (2.1-7)  and  (2.1-9) 
yields  an  expression  for  the  total  free  energy,  F,  of  the 
system  including  the  effect  of  the  centrifugal  field,  i.e. 
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(2.1-10) 

The  condition  for  an  extremum  of  the  free  energy 


can  be  expressed  as 
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(2.1-11) 


The  6  notation  is  used  to  indicate  a  virtual  variation 
of  the  Helmholtz  potential.  The  variation  in  equation 
(2.1-11)  is  subject  to  the  following  constraints: 


T  =  constant 


(2.1-12) 
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dvg  =  vT 


(2.1-13) 
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(2.1-14) 


m  rp 

where  V1  and  represent  the  total  volume  and  total  mole 

numbers  of  the  entire  system.  For  the  closed,  constant- 

T  T 

volume  system  considered,  V  and  in  the  above  equations 

are  constant.  Using  expression  (2.1-10),  the  variation  of 
the  free  energy  in  equation  (2.1-11)  can  be  expressed  as 
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(2.1-15) 


In  general 


6(r^Dm)  =  2r6r  Dm  +  r^SDm 


(2.1-16) 


For  the  liquid-gas  solution  and  the  gas-vapor  mixture,  the 
element  of  mass  Dm  can  be  expressed  as  pDV  or  ^  M . DN . ,  where 

U  1  1 

i=l 

p  is  the  density  and  ML  is  the  molecular  weight  of  component  i. 
Then  equation  (2.1-16)  becomes 

2  n  2 

6(r  Dm)  =  2r6rpDV  +  Z  r  M.  6DN.  (2.1-17) 

•  i  i 

i=l 

Likewise,  for  the  interface,  the  element  of  mass  Dm  can  be 

n  s 

expressed  as  TDA  or  y  M.DN.  .  The  symbol  T.  denotes  the 

L  11  1 

i=l 


28 


adsorption  of  component  i  at  the  interface,  defined  by 


(2.1-18) 


The  total  adsorption  is  the  sum  of  the  individual 
quantities,  i.e. 

n 

r  =  z  r\  (2.1-19) 

i=l 

For  the  interface,  equation  (2.1-16)  becomes 


5(r  Dm)  = 


n  p  c 

2r5r  IDA  +  Z  r  M.  SDN. 

i=i 


(2.1-20) 


Using  expressions  (2.1-15),  (2.1-17),  and  (2.1-20)  in 

equation  (2.1-11)  yields  an  expression  for  the  extremum 
condition  in  terms  of  the  thermodynamic  variables  of  the 
system,  i.e. 


» 

•  /■ 

-  PL  6DVL  + 

2  x  L_T7L  ,  *} 

-  to  r  Srp  DV  +  Z 

•vL 

J 

i — 1 

ii 

•H 

yi 

> 

(y  .  -  —  co  r  M .  )  SDN  . 
l  2  ii 


N.' 

l 


■ 

( 

-  pgsdvg  + 

2  »  G  G  ,  y 
-  u)  rSrp  DV  +  Z 

i=l 

vG 

vG 

/  G  1  2  2,x  rnvr  G 

(y  .  -  to  r  M.  )  SDN . 

l  2  ii 


N. 

l 


+ 


• 

•  * 

n 

ySDA  + 

-  torSr  TDA  +  Z 

i=l 

A  J 

A 

(y . S-  i  w2r2M. ) SDN.S  =  0 

l  2  ii 


N. 

l 


(2.1-21) 


}An  alternate  name  for  the  adsorption  is  the 


superficial  density. 
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The  equilibrium  conditions  are  obtained  from 
equation  (2.1-21)  subject  to  the  constraints  given  by 
equations  (2.1-12),  (2.1-13),  and  (2.1-14).  Since  varia¬ 
tions  in  the  volume  and  surface  area  are  independent  of 
variations  in  the  mole  numbers,  equation  (2.1-21)  can  be 
expressed  as  two  separate  conditions,  i.e. 


» 

* 

-  PL  6DVL  + 

2  x  L  L 

-  03  r  6rp  DV  + 

•  !'vL  ; 

vL 

-  PG  6DV° 


V 


• 

*  f 

+ 

2  p  G  G 

-  03  r6rp  DV  + 

y6DA  + 

« 

i 

vG  J 

A  J 

-  o)  r6r  F  DA  =  0  (2.1-22) 


and 

n 

I 

i— 1 


1  2  2.  L 

—  03  r  )  6DN^  + 


n 

E 

i=l 


/  G  1  2  2.  p  _ .  _  G 

(y .  -  —  03  r  )  6DN . 

l  2  l 


N. 

l 


n 

+  E 
i=l 


,  S  1  2  2.  x  S 

(y .  -  —  03  r  )  6DN . 

i  2  l 


=  0 


(2.1-23) 


The  conditions  of  constraint  given  by  equations  (2.1-13)  and 
(2.1-14)  expressed  in  terms  of  variations  in  volume  and  mole 


numbers  become 
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sdvl  + 


V 


6  DV  =  0 


(2.1-24) 


V 


» 

p 

6DN . L  + 

SDN . G  + 

l 

1 

_  _  L 

G 

N  . 

1 

N  . 

l 

5DN^  =  0,  i=l,2,...,n 


N  . 

l 


(2.1-25) 


Equation  (2.1-22)  represents  the  condition  of 
mechanical  equilibrium.  The  variation  of  the  system  is 
expressed  as  variations  in  the  size  and  position  of  the 
volume  and  surface  elements.  Considering  the  pressure- 
volume  terms  first,  by  simple  differentiation 


» 

P6  DV  =  6 

■ 

PDV  - 

V  J 

v  J 

6PDV 


(2.1-26) 


In  the  equation  above,  6P  represents  the  variation  in 
the  pressure  associated  with  the  change  in  position  of  the 
volume  element  DV.  For  a  volume,  V,  bounded  by  a  surface. 


A, 


PDV  = 


V 


P  5n  DA 


(2.1-27) 


A 
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where  6n  is  the  normal  component  of  the  motion  of  DA, 
positive  when  outward  with  respect  to  V.  Using  equation 
(2.1-27)  in  (2.1-26)  gives 


* 

• 

P  6DV  = 

/\ 

P6n  DA  - 

J  v  J 

A 

6PDV 


(2.1-28) 


In  terms  of  the  liquid-gas  solution,  V^,  and  the  gas-vapor 

Q 

mixture,  V  ,  equation  (2.1-28)  becomes 


• 

» 

PL  6DVL  + 

PG  6DVG  = 

•vL 

vG 

(PG  -  PL) «  6n  DA 
'  A 


6PL  DVL  - 


V 


SPG  dvg 


V 


(2.1-29) 


G  L 

where  (P  -  P  )  is  the  pressure  difference  across  the 

surface  element  DA,  and  6n  is  positive  when  projecting 

G  L 

outward  from  V  into  V  . 

Consider  now  the  surface  terms  in  equation  (2.1-22). 
Again,  differentiating  gives 


■ 

' 

y6DA  =  6 

yDA  - 

A  ■ 

A 

6  yDA 


(2.1-30) 
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where  5y  is  the  variation  in  the  surface  tension  associated 
with  the  component  of  motion  of  the  surface  element  DA  lying 
in  the  surface.  For  the  "dividing  surface".  A,  with 
principal  curvatures  k  and  <  , 


r 


5 


yDA 


s\ 

Y(k^  +  K2)  <$n 


J  A 


A 


(2.1-31) 


where  the  curvatures  are  considered  positive  when  the 

G 

centers  are  inside  V  .  Using  equation  (2.1-31)  in  (2.1-30) 
gives 


» 

• 

y  6DA  = 

y  (k-^  +  k^)  ^n  DA  - 

A  J 

A 

6  yDA 


(2.1-32) 


Finally,  by  means  of  equations (2 . 1-29)  and  (2.1-32), 
the  general  condition  of  mechanical  equilibrium  is  reduced 
to 


( 6 PL  -  pLco2rSr)  DVL  + 


,  G  2  ~  *  t"\t  yO* 

(6P  -  p  10  ror)  DV 


[-(pG-pL)  6n  +  y  (k  +  k9)  6n  -  Too^rfir  -  6y]  DA  =  0 


A 


(2.1-33) 
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Equation  (2.1-33)  is  satisfied  by  the  condition  that  the 
integrands  vanish  throughout  the  system,  i.e. 


dPL  L  2 

dr“  =  p  ^  r 


(2.1-34) 


dPG  G  2 

dF"  =  p  w  r 


(2.1-35) 


-(P 


-  PL) 


A 


Sn  +  y  (k  + 


K2> 


6n  -  T  co  r6r  -  8y  =  0  (2.1-36) 


Separating  the  variations  in  equation  (2.1-36)  into  normal 
and  tangential  motions  of  DA  yields 


(PG  -  PL)a  =  y  (k-l  +  k 2 )  -  r  w2r  ^ 

dn 


(2.1-37) 


and 


dy 

dr 


2 

co  r 


(2.1-38) 


respectively.  If  x  is  defined  as  the  angle  between  the 
normal  and  radial  direction,  i.e. 


dr  =  dn  cos  x 


(2.1-39) 


equation  (2.1-37)  can  be  rewritten  as 


(PG  -  PL)A  =  Y  (  ~  r  ^ r  cos  X  (2.1-40) 


Consider  now  the  conditions  for  chemical  equilibrium 


given  by  equation  (2.1-23).  Since  the  variations  in  mole 
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numbers  are  not  independent  of  each  other,  the  method  of 
Lagrange  undetermined  multipliers  is  used  to  solve  equation 
(2.1-23)  subject  to  the  constraint  given  by  equation 
(2.1-25).  Equation  (2.1-25)  is  multiplied  by  the  constant 
A. ,  and  then  subtracted  from  equation  (2.1-23)  to  give 


/* 


n 

E 

i=l 


(li.L-  ^  co2r2M.-A.)6DN.L  + 
i  2  i  i  i 


n 

E 

i=l 


,  G  1  2  2aa  .  . 

(yi  '  2  “  r  W 


(y.S-  i  u2r2M.  -  A.)  6DN.S  =  0  (2.1-41) 

i  2  i  i  i 

Equation  (2.1-41)  is  satisfied  by  the  condition  that  each 
integrand  vanish,  which  yields 

L  1  2  2„  G  12  2 

y .  -  j  r  Mi  =  *  2  “  r  Mi 

=  y.^  -  =  A.,  i=l ,  2 ,  .  .  .  ,  n  (2.1-42) 

i  2  ii 

The  equilibrium  conditions  for  a  gas-vapor  bubble 
in  a  liquid— gas  solution  subject  to  a  centrifugal  field 
can  be  briefly  summarized  by  the  following  equations: 


n 

+  E 
i=l 
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(2.1-43) 


(pG  “  pL)A  “  y(K2  +  -  ro)2r  cos  x 


(2.1-44) 


dy 


(2.1-45) 


dr 


2  2 
oo  r  M 


i 


2  2 

oo  r  M .  =  A 


i  i 


i=l,  2 


.  .  .  ,  n 


(2.1-46) 


2 . 2  The  Equilibrium  Shape  of  the 
Gas-Vapor  Bubble 

The  equilibrium  shape  of  the  gas-vapor  bubble  at 
the  axis  of  rotation  of  the  system  will  be  deformed  from  a 
spherical  shape  by  the  centrifugal  field.  In  this  section, 
an  expression  for  the  equilibrium  shape  of  the  gas-vapor 
bubble  will  be  derived  from  the  conditions  of  equilibrium 
obtained  in  the  previous  section. 

Since  the  bubble  is  symmetric  about  the  axis  of 
rotation,  the  curve  produced  by  a  plane  section  through 
the  axis  is  used  to  represent  the  shape  of  the  bubble. 
Figure  2.2-1  shows  such  a  plane  section.  The  radial 
direction,  perpendicular  to  the  axis  of  rotation,  is  indi¬ 
cated  by  the  r-axis.  The  r  and  x  intercepts  are  denoted  by 
'a'  and  'b',  respectively.  The  distance  'a'  represents  the 
maximum  radius  of  the  bubble.  The  radius  of  curvature  at 


liquid-gas  solution 
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Figure  2.2-1  A  Plane  Section  Through  the  Axis 

Rotation  of  the  Gas-Vapor  Bubble 
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the  center  of  the  system,  i.e.  at  r=0,  is  denoted  by  'h'. 
The  angle  0  is  that  associated  with  the  slope  of  the  curve. 
The  angle  between  the  normal  to  the  curve  and  the  positive 
x-axis  is  denoted  by  <j) . 

2  7 

Neglecting  the  term  Too  rcosy  in  equation  (2.1-40) 
results  in  Laplace's  equation  of  capillarity,  i.e. 


(P 


G 


Y  (k 2.  +  k 2 ) 


(2.2-1) 


relating  the  pressure  difference  across  the  interface  to 
the  surface  tension  and  principal  curvatures  of  the  inter¬ 
face  at  that  point.  Equation  (2.2-1)  is  further  simplified 
by  assuming  the  surface  tension  to  be  constant  over  the 

g 

interface.  Once  the  pressure  difference  is  expressed  as  a 
function  of  the  distance  r  from  the  axis  of  rotation,  an 
appropriate  expression  for  the  curvature  can  be  substituted 
into  equation  (2.2-1)  to  yield  a  differential  equation  for  the 


7 

The  value  of  the  adsorption,  r,  is  usually  suffici¬ 
ently  small  to  justify  neglecting  the  term  rw2rco  i-n  equation 

(2.1-40).  For  example,  Defay_,  Prigogine,  and  Bellemans  [30] 
obtained  a  value  of  0.62  x  10  1U  kg-mole  m-2  as  the  approximate 
magnitude  of  the  relative  adsorption  of  nitrogen  at  the  plane 
surface  of  a  saturated  solution  of  nitrogen  in  water  with  a 
gaseous  phase  of  nitrogen  and  water  vapor. 

o 

From  the  equilibrium  conditions. 


dy 

dr 


r 


In  view  of  the  small  value  of  the  adsorption,  F,  the  radial 
variation  of  y  can  be  neglected. 
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equilibrium  shape  of  the  bubble. 

Equation  (2.1-34)  gives  the  radial  pressure 
variation  in  the  liquid-gas  solution.  Integrating  equation 
(2.1-34)  from  r=0  where  PL  =  to  a  point  'r'  results  in 

PL  (r)  =  +  -  pLw* 2r2  (2.2-2) 

where  pL  is  constant  due  to  the  assumption  of  liquid 

incompressibility.  For  a  gas-vapor  mixture  of  low  density, 

G  9 

the  pressure  inside  the  bubble,  P  ,  can  be  assumed  constant. 
The  pressure  difference  across  the  interface  at  the  center 
of  the  system  can  be  related  to  the  radius  of  curvature,  ' h' 
at  that  point,  i.e. 


(pG-Pc>A=  It  (2'2-3) 

Using  equations  (2.2-2)  and  (2.2-3),  the  pressure  differ¬ 
ence  across  the  interface  at  any  point 'r' can  be  expressed  as 


.G  L,  2y  1  L  2  2 

(p  - p  >a  =  tr  - 1 p  “  r 


(2.2-4) 


Q 

From  the  equilibrium  conditions, 

2 

=  P  w  r 
dr 

For  a  gas-vapor  mixture  of  low  density,  i.e. 

.G 


pG«  i. 


dP 

dr 


=  0 


.G  . 


implying  that  P  is  approximately  constant, 


' 
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The  shape  of  the  surface  at  any  point  is  character¬ 
ized  by  the  two  principal  curvature,  k.,  and  k^.  For  the 
axi symmetric  surface  being  considered,  is  the  plane 
curvature  of  the  meridional  curve  shown  in  Figure  2.2-1. 

The  second  curvature,  k  ,  acts  in  a  perpendicular  direction 
to  the  x-r  plane.  The  radius  of  curvature  associated  with 
,  denoted  by  in  Figure  2.2-1,  is  equal  to  the  distance 
between  the  curve  and  the  x-axis  along  a  line  normal  to  the 
curve  in  the  x-r  plane.  From  analytical  geometry  [31] ,  the 
sum  of  the  principal  curvatures  can  be  expressed  as 


K1  +  K2  = 


(d£)2  _  r  <^§> 

1  +  ^dx  dx 


r  [1  +  (ff)2] 


3/2 


(2.2-5) 


The  angles  0  and  c p,  positive  in  a  counterclockwise  sense, 
are  related  by  the  equation 


0  =  (p  -  tt/2 


(2.2-6) 


where  0  is  defined  by 


tan  0 


dr 

dx 


(2.2-7) 


Using  equation  (2.2-7)  in  (2.2-6)  gives 


dr 

dx 


cot  (p 


(2.2-8) 


which  upon  differentiating  yields 


' 


d2r 

dx2 


2  d(j) 

CSC  (t)  — - — 

dx 
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(2.2-9) 


Substituting  expressions  (2.2-8)  and  (2.2-9)  into  (2.2-5) 
results  in 


K1  +  K  2 


1  +  cot2(j)  -  (r)csc2(|)  ^ 


r(l  +  cot  0) 


3/2 


(2.2-10) 


which  reduces  to 


K 


1 


+ 


1  d_ 
r  dr 


(rsin4> ) 


(2.2-11) 


Finally,  using  expressions  (2.2-4)  and  (2.2-11)  in 
(2.2-1)  gives  a  differential  equation  for  the  shape  of  the 
gas-vapor  bubble  in  terms  of  the  parameters  'r'  and  4),  i.e. 


1  d_ 
r  dr 


(rsin4>) 


2 

h 


L  2  2 


p  co  r 

2y 


(2.2-12) 


Expression  (2.2-12)  can  be  integrated  to  give 

L  2  3 

,rpcor 

sin  =  —  —  -  (2.2-13) 

h  8Y 

where  the  constant  of  integration  is  zero  since  r=0  at  <J>=0. 

A  relationship  between  the  maximum  radius  of  the  bubble,  'a', 

and  the  radius  of  curvature  at  the  center,  'h'  is  obtained  by 

evaluating  equation  (2.2-13)  at  4)  =  —  where  r=a,  i.e. 

L  2  3 

a  _  ,  ,  p  oj  a 

h  “  1  +  ~8Y  (2.2-14) 


Let  the  shape  parameter  'e'  be  defined  as 
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8  Y 

Then,  from  equation  (2.2-14) 

a  (1  +  e)  h 

and  equation  (2.2-13)  can  be  expressed  as 

3 

sm<J>  =  ( 1  +  e)  R  -  eR 


(2.2-15) 


(2.2-16) 


(2.2-17) 


where  R  is  the  dimensionless  co-ordinate  given  by 


RE-  (2.2-18) 

a 

Equation  (2.2-17)  expresses  the  shape  of  the  bubble  in 
terms  of  R  and  (p ,  using  the  shape  parameter  'e'. 

The  shape  of  the  bubble  is  more  conveniently 
expressed  as  an  explicit  function  for  the  dimensionless 
co-ordinate  X  in  terms  of  R,  where 


X 


X 

a 


From  equation  (2.2-18) 


(2.2-19) 


^  =  -  tan  <j>  (2.2-20) 

dR 

If  tan  <f>  is  expressed  in  terms  of  sin  <j> ,  equation  (2.2-20) 
becomes 


dX 

dR 


-  sine}) 

2  *77 

(1  -  sin  (p) 


(2.2-21) 


Using  expression  (2.2-17)  for  sine})  in  (2.2-1)  results  in 


dX  =  -  R[ (1  +  e)  -  eR  J 

dR 


{1  -  R2  [  (1  +  e)  -  eR2]2} 
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(2.2-22) 


The  solution  to  the  differential  equation  represented  by 
expression  (2.2-22)  can  be  obtained  in  more  than  one  form. 
The  expression  for  dX  can  be  integrated  directly  in  terms 
of  elliptic  integrals  of  the  first  and  second  kind  to  give 
a  general  expression  for  X.  The  result  is  presented  in 
Appendix  A. 

However,  an  alternate  formulation  utilized  by  Lord 
Rayleigh  [32]  in  his  treatment  of  a  rotating  fluid  mass 
subject  to  capillary  tension  expresses  the  solution  as  a 
series  expansion  in  terms  of  the  shape  parameter,  'e'.  Intro¬ 


ducing  the  transformation 


Z  =  1  -  R“ 


into  equation  (2.2-22)  results  in 


(2.2-23) 


dX 

dZ 


1  +  ez 


2 [1  -  (1  -  Z)  (1  +  ez)  ] 


2 ,  v. 


(2.2-24) 


The  solution  to  equation  (2.2-24)  can  be  expressed  as  a 
series  expansion  in  terms  of  powers  of  'e'*  For  the  case  of 
a  small  bubble,  where  e  «  1,  only  the  leading  terms  need 

be  retained. 

An  estimate  of  the  order  of  magnitude  of 'e' for  a 
small  negative  pressure  in  the  solution  at  the  center  of  the 
system  can  be  obtained  using  a  spherical  approximation  for 
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the  bubble  shape.  For  a  given  system,  the  value  of 'e' is 
dependent  upon  the  maximum  radius,  'a',  which  can  be  estimated 
using  the  expression  for  the  critical  radius  of  a  spherical 
vapor  bubble  in  a  pure  liquid  [33],  i.e. 


R 

c 


(2.2-25) 


where  P  is  the  vapor  pressure  of  the  liquid.  For  water 

at  20 °C,  the  surface  tension  is  0.0728  N/m  and  the  vapor 

pressure  is  2.34  kPa  [34].  If  the  liquid  pressure  outside 

the  bubble  is  -100  kpa,  using  equation  (2.2-25),  the  order 

of  magitude  of  'a'  is  1  ym.  For  an  angular  velocity  of  1000 

rad/sec,^  the  order  of  magnitude  of  the  shape  parameter  is 
-9 

10  .  For  larger  negative  pressures,  the  value  of  'e' would 

be  even  smaller.  Thus,  for  negative  pressures  in  the  liquid- 
gas  solution  outside  the  bubble,  the  equilibrium  size  of  the 
gas-vapor  bubble  is  sufficiently  small  to  satisfy  the  condition 
that 

e  «  1  (2.2-26) 

For  such  a  condition,  it  is  mathematically  conven- 

dX 

ient  to  express  the  function  ^  as  a  power  series  expansion 
in  terms  of  'e'.  Equation  (2.2-24)  can  be  rewritten  as 


P 

o 


10 

and  PL 
c 


The  value  of  go  strictly  depends  on  the  values 
for  a  given  system,  i.e. 


co 


P 


V 


of 


. 
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dX 

dZ 


_ 1  +  eZ _ 

2Z  '2  [1  -  2  ( 1-Z )  e  -  Z  ( 1-Z )  e2  ]  '2 


(2.2-27) 


Using  a  binomial  series  [35]  for  the  term 


J  5  [1  -  2 ( 1-Z ) e  -  Z (1-Z) e  ] 


2  -1'2 


(2.2-28) 


in  equation  (2.2-27)  results  in  the  following  expression 

,  _  11 
for  J : 


J  =  1  +  (l-Z)e  +  (|  -  |  Z  +  Z2 ) e2  +  o  (e3) 


(2.2-29) 


Using  expression  (2.2-29)  in  (2.2-27)  results  in 


dX 

dZ 


2Z 


V  [1  +  (l-Z)e  +  (§  -  |  z  +  Z2)e2] 


[e  +  ( 1-Z )  ez ]  +  o  (e  ) 


(2.2-30) 


Equation  (2.2-30)  can  be  rearranged  in  terms  of  powers  of 
Z  to  give 


||  =  |  (1  +  e  +  f  e2)  Z 


"tie2  ZV2  +  0  (e3, 


(2.2-31) 


Transforming  equation  (2.2-31)  back  into  an  expression  in 
terms  of  R  using  equation  (2.2-23)  gives 


11Exoression  (2.2-29)  is  only  strictly  valid  for  the 

n  <  < 

condition  <  -2  (1-Z)  -  Z(l-Z)e  -  1.  Since  0  -  Z  -  1,  this 
condition  is  satisfied  when  e  <<  1. 


dX 

dR 
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=  -  (1  +  e  +  |  e2) 


R ( 1-R2 ) 


3  2  2 

+  j  eZR  (1  -  Rz) 


+  0  (e°) 


(2.2-32) 


Integrating  equation  (2.2-32)  results  in 


X  = 


3  2 


.2.  % 


(1  +  e  +  J  e") ( 1-R")  2 


(1-R2)  '2 


+  O  (e  ) 


(2.2-33) 


where  the  constant  of  integration  is  zero  since  X=0  at  R=l. 
If  only  terms  of  order  'e'  are  retained  in  equation  (2.2-33)  , 
the  curve  produced  by  a  plane  section  through  the  axis  of 
the  bubble  can  be  approximated  by  the  ellipse  given  by 


X 


(1  +  e) 


+  R  —  1 


(2.2-34) 


For  small  bubbles,  where  e  <<  1,  the  bubble  shape  is  essenti¬ 
ally  that  of  the  spheroid  generated  by  rotating  the  curve 
given  by  equation  (2.2-34)  about  the  x-axis,  the  axis  of 
rotation . 

Briefly  summarizing,  this  section  considered  the 
equilibrium  shape  of  a  gas-vapor  bubble  in  a  liquid-gas 
solution  at  the  center  of  a  centrifugal  field.  For  negative 
pressures  in  the  solution,  the  small  size  of  the  bubble  led 
to  the  condition  that  e  «  1,  and  the  curve  produced  by  a 

plane  section  through  the  axis  of  rotation  of  the  bubble  can 
be  approximated  by  the  ellipse  given  by  equation  (2.2-34). 


46 


The  elongation  of  the  bubble  along  the  axis  of  rotation  agrees 
with  intuition.  Due  to  the  centrifugal  field,  the  heavier 
liquid  is  displaced  further  away  from  the  axis  of  rotation, 
causing  the  lighter,  gaseous  phase  to  move  closer  to  the  center 
of  the  system.  Thus,  the  bubble  is  radially  flattened. 

In  this  chapter,  the  equilibrium  state  was  determined 
from  the  extremum  condition  on  the  total  energy  of  the  system. 
The  next  chapter  investigates  the  stability  of  this  equilibrium 
state  by  considering  the  exact  nature  of  the  extremum  condition. 


CHAPTER  3 


THE  NATURE  OF  THE  EQUILIBRIUM  STATE 

The  equilibrium  state  is  characterized  by  a  station¬ 
ary  value  of  the  total  energy  of  the  system.  The  stability 
of  the  equilibrium  state  depends  upon  the  nature  of  this 
extremum.  Minimum  and  maximum  values  of  the  total  energy 
represent  stable  and  unstable  equilibrium  states,  respect¬ 
ively  [35] .  In  order  to  determine  the  nature  of  the  extremum, 
the  total  energy  of  the  system  must  be  examined  for  those 
states  close  to  equilibrium. 

Consider  the  difference  in  the  total  energy  of  the 
system,  AF,  between  the  equilibrium  state,  and  a  new  state 
close  to  but  not  at  equilibrium.  The  subscripts  a  and  3 
will  be  used  to  denote  the  equilibrium  and  non-equilibrium 
states,  respectively.  The  gas-vapor  bubble  in  each  state 
is  characterized  by  a  specific  size  and  shape,  as  shown  in 
Figure  3-1.  From  equation  (2.2-34),  the  size  and  shape  of 
the  equilibrium  bubble  are  dependent  on  the  shape  parameter, 
'e',  and  the  maximum  radius,  'a'.  But  from  equation  (2.2-15)  fe' 
for  a  given  system  is  a  function  of  'a'.  Thus,  'a'  is  the 
obvious  choice  for  the  bubble  parameter  to  which  the  state 
of  the  system  is  referenced.  The  variation  of  the  maximum 
radius  from  the  equilibrium  value  will  indicate  a  displace¬ 
ment  of  the  state  from  the  equilibrium  position.  Let  e  denote 
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Figure  3-1  A  Comparison  of  the  Plane  Sections  of  the 
Gas-Vapor  Bubble  in  the  Equilibrium  and 
Non-Equilibrium  States 
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the  difference  in  the  maximum  radius  of  the  bubble  between 
the  equilibrium  and  non-equilibrium  states,  i.e. 

e  i  ag  -  aa  (3-1) 

where  e  can  be  positive  or  negative.  Using  equation  (3-1) 
in  expression  (2.2-15)  for  ' e'  allows  the  value  of  the 
shape  parameter  in  both  states  to  be  related,  i.e. 


/ .  £ »  3 

en  =  e  (1  +  — ) 

3  a  a 


(3-2) 


a 


In  chapter  2  an  expression  was  developed  for  the 
total  energy  of  the  equilibrium  system,  i.e. 
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(3-3) 
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For  the  new  state"*-  slightly  displaced  from  equilibrium  the 
total  energy  of  the  system  can  likewise  be  expressed  as 


F  = 
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PB  DvL  +  ih 
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(3-4) 


Subtracting  equation  (3-3)  from  (3-4)  gives  an  expression 
for  the  difference  in  the  total  energy  of  the  system  between 
the  new  and  equilibrium  states,  i.e. 


Although  the  system  in  the  displaced  state  is  not 
in  equilibrium  as  a  whole,  the  individual  phases  can  be 
assumed  to  exist  in  a  state  of  "local  equilibrium"  for  which 
the  equilibrium  function  for  the  Helmholtz  potential,  F,  is 
valid  [ 37 ]  . 
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(3-5) 
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Consider  first  the  pressure-volume  terms  in  equation 
(3-5)  ,  i.e. 


» 

• 

* 

-  PeL  dvl  - 

-  p£  DVL  + 

-  PB  DvG  - 

V3  ' 

vL 

a 

^  J 

(3-6) 

Introducing  the  change  in  volume  AV,  defined  by 


AV  E  -  VG  =  VL  - 

3  oi  a  3 


(3-7) 


into  expression  (3-6)  results  in 
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(3-8) 


Since  the  bubble  is  very  small  compared  to  the  surrounding 
liquid-gas  solution,  the  pressure  distribution  in  the  liquid- 
gas  solution  in  both  states  can  be  assumed  equal,  i.e. 


(3-9) 
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Condition  (3-9)  reduces  expression  (3-8)  to 


-  PG)  DV°  + 
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-  PL)  DV 
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AV 


(3-10) 


Now  consider  the  chemical  and  field  potential  terms 
in  equation  (3-5)  ,  expressed  by 
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(3-11) 


Since  the  number  of  molecules  associated  with  the  bubble 
interface  is  much  less  than  the  number  of  molecules  in  the 
gas— vapor  mixture  and  the  liquid— gas  solution r  the  last  two 
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terms  can  be  neglected  in  expression  (3-11)  as  small  compared 
to  the  other  terms.  From  the  equilibrium  conditions, 

UiL  -  T  “rM  =  y.G  -  ~  u2r2M.  =  A  ,  i=l,2, . . . ,n  (3-12) 

a  1  ^-a  2  11 

where  A^  is  a  constant.  Using  equation  (3-12)  expression 
(3-11)  can  be  written  in  the  form 
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(3-13) 


Since  the  temperature  and  pressure  distribution  in  the 
liquid-gas  solution  in  both  states  are  equal,  it  follows 
that 

U.L  =  ,  i=l , 2 , . . . ,n  (3-14) 

16  a 

From  equations  (3-12)  and  (3-14)  it  is  evident  that  the 
first  term  in  expression  (3-13)  is  equal  to  zero.  Thus, 
the  contribution  of  the  chemical  and  field  potential  terms 
in  the  expression  for  AF  reduces  to 
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(3-15) 


Finally ,  the  surface  terms  in  equation  (3-5) 
consist  of 


Y0  DA  - 
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Y  DA 
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(3-16) 


Since  the  surface  tension  is  assumed  constant,  expression 
(3-16)  reduces  to 


y  AA 


(3-17) 


where 


AA  =  AD  -  A 

3  a 


(3-18) 


Substituting  the  reduced  expressions  for  (3-6)  , 
(3-11)  ,  and  (3-16)  back  into  equation  (3-5)  for  AF  results 
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The  pressure  and  the  chemical  potentials  in  the 
gas-vapor  mixture  in  the  equilibrium  state  are  functions 

Q 

of  the  maximum  radius,  'a',  and  the  mole  numbers,  N  ,  inside 
the  bubble.  The  new  state  is  considered  sufficiently  close 
to  the  equilibrium  state  to  enable  the  pressure  and  chemical 
potentials  in  the  gas-vapor  mixture  in  the  new  state  to  be 
expressed  by  Taylor  series  expansions  around  the  equilibrium 
values, ^  i.e. 
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'The  notation 
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G  G 

indicates  that  the  partial  derivative  of  the  pressure  P  (a,N  ) 
with  respect  to  'a'  is  evaluated  for  the  equilibrium  values 


a 


and 


N 


>06 
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From  the  equilibrium  conditions,  the  leading  term  in 
equation  (3-21)  can  be  expressed  as 
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(3-22) 


If  the  higher  order  terms  are  neglected  in  equations  (3-20) 
and  (3-21) ,  and  the  resultant  expressions  substituted  into 
equation  (3-19)  ,  the  expression  for  AF  becomes 


AF  = 


-  <pr  - 


av 


a 


P  )  DV  +  yAA 

a 


+ 


V 


a 


3P 

3a 


(e) 


3P 


a  ,NG 

a  ~a 


n 

E 

i=l  3N/ 
i 


(N, 


a  / 

a  ~a 


G 


«iG> 

a 


DV 


+ 


r  ,  G 

?  3PG 

3P 

( £  )  +  E  r 

3  a 

i=l  3N. 

1 

AV 

a  ,NG 

a  ~a 

(Ni  - 

3 

G 


^  ,  N 

a  ~a 


NiG, 

a 


DV 


n 

+  E 
i=l 


r  3y. 


3a 


N  . 


X3 


n  3y  • 

(e)  +  E 

j=l  3N. 

G  3 

a 

a  ~a 


(N  G 

]6 


n.g) 

J  a 


a  ,NG 

a  ~a 


DN. 

l 


(3-23) 
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G  G 

Since  the  quantities  e  and  (N«  -  N  )  can  be  made  arbitrarily 

~  p)  ~  uc 

small,  only  the  first  two  terms  will  be  retained  as  signifi¬ 
cant  in  equation  (3-23) .  Then  the  expression  for  AF  reduces 
to 


AF  = 


-  (P  - 
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P  )  DV  +  yAA 
a  1 


(3-24) 


AV 


Using  expression  (2.2-4)  for  the  equilibrium  pressure 
difference  across  the  interface  in  equation  (3-24)  results  in 
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Equation  (3-25)  can  be  separated  into  three  terms,  i.e. 


AF  =  -J1  +  J2  +  J3 
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The  evaluation  of  ,  J^,  and  requires  the  use  of 
analytical  expressions  for  the  area,  volume,  and  a  term 
similar  to  the  moment  of  inertia  of  the  bubble  in  both  states. 
The  moment  of  inertia,  I,  of  an  axi-symmetric  body  about  the 
cylindrical  axis  is  defined  by 


pr2 * 4  DV 


V 


(3-30) 


The  particular  moment  of  inertia  used  to  evaluate  expression 
(3-28)  for  J 2  is  that  of  a  spheroid  of  liquid  identical  in 
size  and  shape  to  the  gas-vapor  bubble,  and  of  constant 
density  pL.  In  the  previous  chapter  it  was  shown  that  for  a 
small  bubble  at  negative  pressures  the  condition  e  «  1 

enabled  the  equilibrium  bubble  shape  to  be  approximated  by 
the  spheroid  represented  by 


- x  -T  +  =  1  (3-31) 

(l+e)2a  a2 

2  3 

where  terms  of  power  e  and  higher  have  been  neglected. 

For  this  approximation,  the  area,  volume,  and  moment  of 

inertia  required  above  are  given  by  the  following  expres- 
4 

sions ; 


2Since  the  original  expression  for  the  shape  of  the 

bubble  has  neglected  terms  of  power  e2  and  higher,  these  terms 
will  not  be  retained  in  the  following  development. 

4The  derivations  of  the  expressions  for  the  area, 
volume,  and  moment  of  inertia  are  presented  in  Appendix  B. 
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(3-34) 


The  expressions  above  are  in  a  strict  sense  only  valid 

for  the  equilibrium  state.  However,  the  new  state  will  be 

considered  sufficiently  close  to  equilibrium  to  enable  the 

equilibrium  expressions  for  the  shape,  area,  volume,  and 

moment  of  inertia  to  be  used  in  the  displaced  state. 

Now  consider  the  term  J.,  .  Since  both  y  and  h  are 

1  a 

constant  for  a  given  system,  the  integral  in  equation  (3-27) 
reduces  to 


(3-35) 


Using  expression  (3-33) ,  the  difference  in  volume  can  be 
expressed  as 
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Substituting  expressions  (3-1)  and  (3-2)  for  a^  and  e^  in 
equation  (3-36)  ,  and  collecting  terms  results  in 
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Equation  (2.2-16)  expressed  the  equilibrium  radius  of 

curvature  at  the  center  in  terms  of  the  maximum  radius 

and  shape  parameter.  Substituting  equations  (2.2-16)  and 

2 

(3-37)  back  into  (3-35),  and  neglecting  terms  of  power  e  , 
yields  an  expression  for  in  terms  of  the  equilibrium 
parameters  a^  and  e  ,  and  the  radial  displacement  e,  i.e. 

J,  =  8Trya  (1  +  3e  )  e  +  8tty  (1  +  6e  )  e2 

1  1  a  a  a 

+  ^  (1  +  21e  )  —  +  0 ( e4 )  (3-38) 

j  a  aa 

Next,  consider  the  term  J^.  Using  expression  (3-30) 
in  equation  (3-28)  gives 


J 


2 


(3-39) 


where  the  density  of  the  spheroid  considered  is  that  of 
the  liquid,  pL.  Using  expression  (3-34)  for  the  moment 
of  inertia  terms  in  (3-39)  results  in 
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(3-40) 


Again  substituting  expressions  (3-1)  and  (3-2)  for  a^  and 

e0  in  equation  (3-40) ,  and  collecting  terms  gives 

3 
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Using  the  substitution 
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equation  (3-41)  becomes 
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-  [  (5  +  8ea)  aa  e  +  2(5  +  14ea)e2 


+  2(5  +  2 8e  )  |-  ]  +  0(e4) 

DC  ct 

a 


(3-43) 


Neglecting  powers  of  e  reduces  expression  (3-43)  for  to 


32Trye  ?  3  . 

J2  =  +  2,  +  2  §-)  +  0(c) 

a 


(3-44) 


Finally,  consider  the  term  J^.  Using  expression 
(3-32)  for  the  area  in  (3-29)  results  in 


J3  =  4 tty  [aB2  (1  +  f  eB)  -  aa2  (1  +  f  ej  ]  (3-45) 

Substituting  expressions  (3-1)  and  (3-2)  for  a0  and  e0  in 

p  p 

the  above  equation  yields  the  following  expression  for  J^: 
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J3  =  87rYaa  (1  +  I  ea>  £  +  4,IY  (1  +  F  ea)e2 


+  r-  7TYea  |-  +  0(e4) 

a 


(3-46) 


Returning  to  expression  (3-26),  and  substituting 
for  J^,  J2 ,  and  from  equations  (3-38),  (3-44),  and 

(3-46)  ,  respectively,  results  in  an  expression  for  AF 
close  to  equilibrium,  i.e. 


AF  =  4wyaa2  [  - (|— )  2  -  §  (1  +  3ea)  (|— )  3]  +  0(e4) 

a  a 


(3-47) 


where  ( — )  «  1. 
a 

a 


Expression  (3-47)  retains  terms  of  0(e  ).  As  a 
means  of  determining  how  many  terms  must  be  retained  in 
equation  (3-47)  ,  the  expression  for  AF  when  no  field  is 
present  can  be  compared  to  that  obtained  by  Ward  et  ad.  [38] 
for  the  case  of  positive  pressures  in  the  liquid-gas  solution 
When  the  angular  velocity  is  zero,  the  shape  parameter  is 
zero,  and  equation  (3-47)  reduces  to 


2  p  2  2  f  3  4 

AF  =  47rya  ^  [  -  (f-)  -  f  (§-)  ]  +  0(e  ) 

a  aa  J  aa 


(3-48) 
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Ward  et  al.  obtained  the  following  expression3  for  the 
difference  in  the  free  energy  of  the  equilibrium  and  non¬ 
equilibrium  systems: 


AF  =  4iry 


2  rA 

3  R  ’ 

c 


(3-49) 


where  Rc  and  R  are,  respectively,  the  radii  of  the 
spherical  gas-vapor  bubble  in  the  equilibrium  and  non¬ 
equilibrium  states.  For  means  of  comparison,  the  expression 


R  =  R  +  e 
c 


(3-50) 


can  be  substituted  into  equation  (3-49)  to  give 


AF 


4TryR 


2 

3 


(— )  3] 

V 


(3-51) 


Comparing  equations  (3-48)  and  (3-51)  ,  and  noting  that 
for  zero  angular  velocity,  a  =  R  ,  it  is  evident  that 

0(r  O 

5 

Ward  et  ad.  obtained  the  following  expression  for 
the  change  in  the  Helmholtz  potential,  AF  ,  due  to  the  forma¬ 
tion  of  a  spherical  gas-vapor  bubble  of  °  radius  R  close  to 

the  equilibrium  radius,  R  ,  in  a  liquid-gas  solution: 

c 

AFq  =  4 TTy  (R2  -  |  |-) 

c 

In  order  to  obtain  the  difference  in  the  Helmholtz  potential 
between  the  equilibrium  and  non-equilibrium  states,  the 
reversible  work,  WR,  for  the  formation  of  the  equilibrium 
bubble  must  be  subtracted  from  the  equation  above.  The 

reversible  work  is  obtained  from  the  equation  above  evaluated 
at  R  =  R  ,  i . e . 

c  2 

4ttyR 

WR  = 


3 
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expression  (3-47)  for  AF  reduces  to  that  of  Ward  et  al.  for 
the  no  field  condition.  Terms  of  O(e^)  can  be  neglected  in 
equation  (3-47) ,  and  the  expression  for  AF  reduces  to 

AF  =  47TYaa2  [  -  (|— )  2  -  |  (1  +  3ea)  (|— )  3]  (3-52) 

a  a 

Since  equation  (3-52)  represents  the  change  in 
the  free  energy  of  the  system  as  a  function  of  the  radial 
displacement  e ,  it  can  be  used  to  determine  the  nature  of 
the  extremum  at  equilibrium,  i.e.  at  e=0.  Differentiating 
equation  (3-52)  twice  with  respect  to  e  gives 


d  (Af) 
de 


4TTYa 

a 


2(1  +  3ea)(f-)2] 

a 


(3-53) 


and 


d  (AF)  . 

- ^ - -  =  4  tty 

de 


[  -  2  -  4  (1  + 


3e  )  (  — 
a  a 


-)  ] 


a 


(3-54) 


d  ( AF ) 

Setting  — ^ equal  to  zero  yields  the  stationary  value,  e=0, 

2 

for  which  d  (Af)  is  less  than  zero.  Thus,  the  extremum 

de2 

for  the  free  energy  of  the  system  at  equilibrium  is  a 

maximum.  Figure  3-2  shows  the  graph  of  equation  (3-52)  in 

_  6 

the  region  close  to  e=0  for  e  =  10 

Since  the  free  energy  of  the  system  is  a  maximum 
for  e=0,  the  equilibrium  state  of  the  gas-vapor  bubble  is 
unstable.  The  condition  of  instability  leads  to  the  concept 
of  a  critical  size  for  the  bubble,  characterized  by  the 
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AF 


Figure  3-2  The  Behavior  of  the  Helmholtz  Potential 

in  the  Neighborhood  of  the  Equilibrium 
State  for  Which  e  Equals  Zero.  The  value 
of  the  shape  parameter  'e'  is  1  iq-6 
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equilibrium  value  of  the  maximum  radius,  a  .  If  a  bubble 

a 

is  created  of  less  than  equilibrium  size,  the  decrease  in 
energy  accompanying  bubble  contraction  causes  the  bubble 
to  spontaneously  collapse.  Likewise,  if  a  bubble  is  created 
of  greater  than  equilibrium  size,  the  decrease  in  energy 
accompanying  bubble  expansion  causes  the  bubble  to  spontan¬ 
eously  grow.  In  this  sense,  the  equilibrium  size  represents 
a  critical  value  which  must  be  exceeded  if  the  bubble  is  to 
grow . 


Stability  considerations  have  shown  that  the  free 


energy  of  the  gas-vapor  bubble  is  a  maximum  for  the  equili¬ 
brium  state.  The  reversible  work  of  formation,  W  ,  for  the 

I\ 

equilibrium  gas-vapor  bubble  represents  the  increase  in  the 
free  energy  of  the  equilibrium  system  due  to  the  formation 
of  the  gas-vapor  bubble.  An  expression  can  be  obtained  for 
W  by  comparing  the  Helmholtz  potential,  F  ,  of  the  equili- 
brium  system  with  the  Helmholtz  potential  of  the  liquid-gas 
solution  without  the  bubble  present.  The  difference  in  free 
energy  between  the  two  systems  represents  the  work  required 
to  form  the  equilibrium  gas-vapor  bubble.  If  the  subscript  v 
is  used  to  denote  the  equilibrium  system  consisting  of  only 
the  liquid-gas  solution,  the  expression  for  F^  is  expressed 


by 


F 


v 


dvl  +  1 

i=l 


,  L  1  2  2  v  L 

(yi  -  j  w  r  Mi)  DNi 


L 


v 


(3-55) 
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The  condition  of  chemical  equilibrium  in  the  liquid-gas 
solution  requires  that 

L  12  2 

^  r  M.  =  A^  ,  i=l,2,...,n  (3-56) 

v  1  1v 


where  A^  is  a  constant.  Using  expression  (3-56)  in 

v 

(3-55)  gives 


F 

v 


n 

E 

i=l 


(3-57) 


Equation  (2.1-42)  for  the  condition  of  chemical  equilibrium 
in  the  system  with  the  bubble  present  can  be  used  to  reduce 
expression  (3-3)  for  to 


F 

CL 


* 

• 

_  L  L 

G  ttG 

1 

•x) 

1 

o 

< 

1 

+ 

-  P  DV  + 

a 

a 

L 

V 

V 

a, 

CL 

yDA  + 


A 


n 

E 

i=l 


A  .N. 


l  l 


(3-58) 


Subtracting  expression  (3-57)  for  F  from  (3-58)  for  F  gives 

v  AX 

the  following  expression  for  the  reversible  work  of  formation 
for  the  equilibrium  gas-vapor  bubble: 
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WR  = 

* 

-  dvl  + 

•  f 

-  PG  dvg  - 

a 

■ 

vL 

vG 

a 

a 

n 


yDA  +  £ 


[A.  N. 

•  ~i  1  1 

i=l  a  a 


A  .  N  .  ] 


A 


Noting  that 


vL  +  VG 

a  a 


and 


N  . 

l 


T 

a 


(3-59) 


(3-60) 


(3-61) 


in  equation  (3-59)  results  in 


-  (K,  - 


a 


Py)  dvl  + 


V 


a 


-  (P^  - 


a 


PL)  DVG 

v 


V 


a 


+ 


yda 


n 

+  Z  ( A  .  —  A  ,  ) 
i=l  1a  1v 


A 


(3-62) 


Since  the  bubble  is  very  small  in  comparison  to  the 


surrounding  solution,  the  pressure  distribution  in  the 


liquid-gas  solution  in  both  systems  can  be  assumed  equal. 


70 


i .  e . 


P 


L 

a 


(3-63) 


The  equality  of  the  pressure  distribution  and  temperature 
in  the  liquid-gas  solution  in  both  systems  implies  that 


,  1  1,2,...,  n 


(3-64) 


Using  expressions  (3-12)  and  (3-56)  in  equation  (3-64) 
results  in 


A^  =  A^  ,  i=l,2,...,n  (3-65) 

a  v 

Equations  (3-63)  and  (3-65)  are  used  to  reduce  expression 

(3-62)  for  W  to 
K 


* 

-  (pg  -  pL)  dvg  + 

a  a 


yDA 


J  A 


(3-66) 


Using  expression  (2.2-4)  for  the  pressure  difference 
across  the  equilibrium  interface,  and  noting  that  y  is 
constant  in  expression  (3-66)  results  in 


1 

2 


L  2  2 

P  co  r 


) DV  +  yA 


(3-67) 
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or  upon  rearranging. 


2yV 

h 

ot 


L  2 


P  r 


DV  +  yA 


V 


(3-68) 


where  V  and  A  refer  to  the  equilibrium  bubble.  Substitu¬ 
ting  expression  (3-33)  for  V,  and  expression  (2.2-16)  for 
h^  in  the  first  term  of  equation  (3-68)  gives 


2yV 
h 


8tty  a 


a 


a 


(1  +  2e  ) 

a 


(3-69) 


where  powers  of  ea  have  been  neglected.  If  the  integral  in 
the  second  term  of  equation  (3-68)  is  evaluated  using  expres¬ 
sion  (3-34)  for  the  moment  of  inertia,  the  resultant 
expression  is 


1 

2 


pLr 2  DV 


V 


47TpLu)2a  5  (1  +  e) 

a 

15 


(3-70) 


Using  expression  (3-42)  in  (3-70)  ,  and  neglecting  powers  of 

e2  results  in 

a 


1 

2 


2 

0) 


pLr2  DV 


32Tiya  2e 

a  a 

15 


(3-71) 


V 
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Substituting  expression  (3-32)  for  A  in  the  last  term  of 
equation  (3-68)  gives 

yA  =  47rya  2  (1  +  \  e  )  (3-72) 

a  3  a 

Finally,  using  expressions  (3-69)  ,  (3-71)  ,  and  (3-72)  in 

(3-68)  yields  an  expression  for  the  reversible  work  of 
formation  for  the  equilibrium  gas-vapor  bubble,  i.e. 

/i  2 

4TTya 

WR  =  -y-  (1  -  |  ea)  (3-73) 

In  summary,  the  stability  analysis  determined  that 
the  free  energy  of  a  gas-vapor  bubble  in  a  liquid-gas  solution 
at  negative  pressure  is  a  maximum  for  the  equilibrium  state. 
Consequently,  the  equilibrium  state  is  unstable,  and  the 
equilibrium  size  of  the  bubble  represents  a  critical  value. 

A  gas-vapor  bubble  must  exceed  the  critical  size  in  order  to 
spontaneously  grow.  The  reversible  work  required  to  form  a 
gas-vapor  bubble  of  critical  size  is  expressed  by  equation 
(3-73) . 

In  the  next  chapter,  an  expression  is  developed  for 
the  critical  size  of  a  gas-vapor  bubble  in  a  dilute,  liquid- 
gas  solution  at  negative  pressure  in  terms  of  the  properties 


of  the  solution. 


CHAPTER  4 


THE  CRITICAL  SIZE  OF  A  GAS-VAPOR  BUBBLE  IN  A 
DILUTE  SOLUTION  AT  NEGATIVE  PRESSURES 

In  chapter  2  an  expression  for  the  equilibrium 
shape  of  the  gas-vapor  bubble  was  obtained  from  the 
conditions  of  equilibrium.  For  the  case  of  negative 
pressures  in  the  solution  at  the  center  of  the  system,  e  <<  1, 
and  the  bubble  was  approximated  by  the  spheroid  represented 
by  equation  (2.2-34).  In  the  case  of  a  spherical  bubble 
when  no  field  is  present,  the  size  is  expressed  in  terms  of 
the  equilibrium  or  critical  radius  [39] .  For  a  spheroidal 
bubble,  equation  (2.2-34)  expresses  the  equilibrium  size  and 
shape  as  a  function  of  the  maximum  radius  'a'. 

Since  the  bubble  size  and  shape  is  an  equilibrium 
property  of  the  system,  it  is  appropriate  to  develop  an 
expression  for  'a'  in  terms  of  the  properties  of  the  liquid- 
gas  solution.  However,  the  value  of  'a'  is  implicitly  expressed 
by  the  cubic  equation  given  by  (2.2-14),  which  depends  upon 
the  radius  of  curvature  at  the  center  of  the  bubble,  'h'. 

Thus,  '  h'  is  the  basic  parameter,  which  when  evaluated  in 
terms  of  the  properties  of  the  solution,  enables  the  size  and 
shape  of  the  equilibrium  gas-vapor  bubble  to  be  determined. 

As  the  bubble — and  hence  the  shape  parameter — becomes  smaller 
and  smaller,  the  deviation  from  sphericity  becomes  practically 
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negligible.  For  a  perfectly  spherical  bubble,  the  radius 
and  radius  of  curvature  are  equal.  In  practice, 
since  e  <<  1,  evaluation  of  '  h'  in  terms  of  the  solution  prop¬ 
erties  will  be  sufficient  to  describe  the  equilibrium  size 
of  the  bubble. 

Up  until  now,  the  analysis  has  been  of  a  general 
nature,  i.e.  the  equilibrium  and  stability  considered  the 
case  of  a  multicomponent  liquid-gas  solution.  At  this  point, 
the  development  will  be  restricted  to  the  specific  case  of 
a  dilute  solution  of  a  single  gas  in  a  liquid  solvent.  The 
term  dilute  implies  that  the  number  of  solute  molecules  is 
much  less  than  the  number  of  solvent  molecules,  and  conse¬ 
quently,  any  interaction  among  the  solute  molecules  may  be 
ignored  [40] .  The  concentration,  C,  of  the  solution  is 
defined  as  the  ratio  of  the  number  of  moles  of  solute  to  the 
number  of  moles  of  solvent,  i.e. 


where  the  subscripts  1  and  2  refer  to  the  liquid  and  gas 
components,  respectively.  For  a  dilute  solution, 

C  <<  1  (4-2) 

In  addition,  it  is  assumed  that  the  gas  and  vapor  inside  the 
bubble  behave  as  an  ideal  gas  mixture.  Using  appropriate 
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expressions  for  the  properties  of  a  dilute  solution  and  an 

ideal  gas  mixture,  an  expression  for  ' h'  can  be  developed  in 

terms  of  the  liquid  phase  properties. 

The  value  of  ' h'  is  obtained  from  Laplace's  equation 

of  capillarity,  given  by  (2.2-1),  evaluated  at  the  center 

2y 

of  the  system  where  the  bubble  curvature  is  ^  .  Rearranging 
equation  (2.2-3)  gives 


h  = 


2y 


<pG  -  pl> 

c  c 


(4-3) 


Integrating  equation  (2.1-34)  from  r=£  where  P  =  Pq  to  r=0 

yields  an  expression  for  PL,  the  pressure  in  the  liquid  at 

c 


the  center  of  the  system,  i.e. 


PL  =  P  -  i  pL(A2  (4-4) 

c  o  2 

Q 

It  remains  to  obtain  an  expression  for  P  in  terms  of  the 

c 

properties  of  the  dilute  solution. 

The  partial  pressures  of  the  gas  and  vapor  components 
inside  the  bubble  can  be  expressed  in  terms  of  their  respect¬ 
ive  mole  fractions,  and  X i.e. 

PG  E  PGX.  ,  i=l , 2  (4-5) 

1  1 

where 

N  . 

l 


X  . 

l 


i=l ,  2 


(4-6) 
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and 

N  E  N1  +  N2  (4-7) 

The  total  pressure,  P^,  is  the  sum  of  the  partial  pressures, 

i .  e . 

PG  =  PG  +  PG  (4-8) 

The  expressions  for  the  partial  pressures  are  obtained 
from  the  condition  of  chemical  equilibrium  across  the 
interface,  i.e. 


1 

2 


2  2 
u)  r  M 


1 

2 


2  2 

oo  r  , 


i=l ,  2 


(4-9) 


The  general  expression  for  the  chemical  potential  of  each 
component  in  an  ideal  gas  mixture  is 


Mi (P /T) 


(Pr,T) 


+  kT  in 


(4-10) 


G 

where  1  k'  is  Boltzmann's  constant  and  y  .  (P  ,T)  is  the 

Oi  r 

chemical  potential  of  the  pure  component  at  seme  reference  pressure  P^_  [41]  . 

For  a  dilute  solution  of  concentration  C,  the 
chemical  potentials  of  the  solvent  and  solute  can  be 
expressed  respectively  as 


and 


]i±  (P,T)  =  yQ^  (P,T)  -  kTC 
y ^  (P  ,T)  =  \p  (P  ,T)  +  kT  in C 


(4-11) 


(4-12) 
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where  Pq^(P,T)  is  the  chemical  potential  of  the  pure  liquid 
component,  and  i|j(P,T)  is  a  function  of  P  and  T  [42].  For 
an  incompressible  solvent  with  specific  volume  vL,  a 
reference  pressure,  P^_,  can  be  introduced  into  expression 
(4-11)  to  give 


P1(P,T)  =  (Pr,T)  +  VL  (P  -  Pr)  -  kTC 


(4-13) 


Now  evaluating  equation  (4-9)  at  the  interface 
yields  the  condition  of  equality  of  chemical  potentials 
across  the  interface,  i.e. 


vT1  =  ,  i=l,2 


(4-14) 


Substituting  the  appropriate  expressions  for  the  chemical 
potentials  in  the  equation  above  results  in  expressions 
for  the  partial  pressures  in  terms  of  the  properties  of  the 
dilute  solution. 

Consider  first  component  1.  Choosing  P  ,  the  vapor 
pressure  of  the  liquid,  as  the  reference  pressure,  and 
using  expression  (4-10)  and  (4-13)  in  (4-14)  gives 


U01  (pv'T>  +  vL  <pL  -  V  -  kTC 


P^X 

=  y  ®  (P  ,T)  +  kT  ln 

V 


(4-15) 
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Noting  that  y^  (Pv#T)  =  y  ^  (Pv,T)^  and  rearranging  equation 
(4-15)  results  in  an  expression  for  the  partial  pressure  of 
the  vapor  inside  the  bubble,  i.e. 

G  L 

p  xi  =  pv  exP  (p  -  Pv)  -  C]  (4-16) 

Now  consider  component  2.  Using  expressions  (4-10) 
and  (4-12)  in  (4-14)  gives 

p^x 

kT£n  (-^-^)  =  ip  (PL  ,T)  -  yQ2  (Pr,T)  +  kT  Z  n  C 

(4-17) 

Ij  g 

The  functions  ip  (P  ,T)  and  ^02^Pr'T^  are  unknown  in  terms 

of  the  properties  of  the  dilute  solution.  Therefore,  a 

reference  condition  will  be  used  to  eliminate  the  term 
L  G 

^(P  ,T)  -  Po2^Pr'T)  from  expression  (4-17).  Since  the 
equilibrium  size  and  shape  of  the  bubble  pertain  to  negative 
pressures  in  the  surrounding  solution,  the  reference  condi¬ 
tion  must  remain  valid  when  the  liquid  near  the  center  of 
the  system  enters  a  state  of  tension. 

Ward  et  al.  [44]  considered  a  reference  system 
consisting  of  a  liquid  saturated  with  the  gas  component 
across  a  flat  interface.  The  pure  gas  and  liquid  phases 

■*"The  vapor  pressure  is  the  pressure  of  the  vapor 
in  equilibrium  with  a  liquid  phase  for  a  flat  interface 
at  a  given  temperature  [43]  .  Since  the  system  is  in 
equilibrium,  it  follows  that  the  chemical  potentials  in 
each  phase  are  equal. 
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were  separated  by  a  non-rigid  membrane  permeable  only  to 
the  gas  component.  For  such  a  system  at  equilibrium,  the 
equality  of  the  chemical  potential  of  the  gas  component 
across  the  interface  enabled  the  following  substitution  to 
be  made  in  equation  (4-17)  with  the  reference  pressure 
chosen  to  be  PL : 

iMPL,T)  -  y 02  (PL,T)  =  -kT  £n  Cq  (4-18) 

where  Cq  is  the  saturation  concentration  of  the  gas  in  the 
liquid  phase.  However,  such  a  reference  condition  is  not 
appropriate  to  a  system  where  the  dilute  solution  outside 
the  bubble  is  subject  to  negative  pressures.  Since  the 
membrane  in  the  reference  system  is  flat,  a  state  of  tension 
in  the  solution  would  require  the  pressure  in  the  gas  phase 
to  also  be  negative,  which  is  impossible.  Thus,  the  refer¬ 
ence  concentration  used  by  Ward  et  al_.  in  their  treatment 
of  an  equilibrium  gas-vapor  bubble  is  not  defined  for  the 
case  of  negative  pressures  in  the  solution  outside  the  bubble. 

Consider  instead  the  reference  system  given  in 
Figure  4-1,  where  the  non-rigid  membrane  at  r=£  is  only 
permeable  to  the  gas  component.  The  liquid  and  gas  sub¬ 
systems  are  denoted  by  the  superscripts  LR  and  GR,  respect¬ 
ively.  Since  the  geometry,  angular  velocity,  and  reference 
pressure  in  the  primary  and  reference  systems  are  the  same, 
the  pressure  distributions  in  the  liquid-gas  solution  are 
identical.  The  condition  of  chemical  equilibrium  for  the 
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Primary 

System 
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System 


Figure  4-1  A  Comparison  of  the  Primary  and  Reference 

Systems 
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gas  component  in  the  reference  system  requires  that 

LR  1  2  2  GR  1  2  2 

U2  “  2  W  r  M2  =  y2  ”  2  °°  r  M2  (4-19) 

Specifically,  consider  condition  (4-19)  for  a  point  *r' 
in  the  liquid  solution  and  the  point  r=&  in  the  gas  at  the 
membrane.  Using  the  expressions  for  the  chemical  potential 
of  the  gas  component  in  each  subsystem  in  equation  (4-19) 
gives 


<MPLR,T)  +  kTSLnC  -  i  U2r2M,  =  uJt  (P  ,T)  -  i  u>2«.V 

s  Z  Z  02  o  2  2 

(4-20) 

where  P  is  the  reference  pressure  at  the  membrane,  and  C 
o  s 

is  the  equilibrium  concentration  at  any  point  '  r'  in  the 

reference  system.  Rearranging  equation  (4-20)  and  noting 
LR  L 

that  P  is  equal  to  P  ,  results  in 

1MPL,T)  -  U02(Po-T)  =  kTJln  (i— )  +  |  to2  (r2  -  £2)M2  (4-21) 


Knowledge  of  the  equilibrium  distribution  of  the 

gas  in  the  reference  system  allows  the  concentration 

at  ' r'  to  be  related  to  the  concentration  C  at  the  membrane. 

so 

The  expression  for  the  equilibrium  concentration,  C  ,  is 
obtained  from  the  conditions  of  equilibrium  in  the  dilute 
solution  in  the  reference  system,  i.e. 


dP 


LR 


dr 


L  2 

p  co  r 


(4-22) 
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and 


LR  1  2  _  ,R 

y2  2WrM2  ^2 


(4-23) 


p 

where  is  a  constant.  Using  the  expression  for  the 
chemical  potential  of  the  solute  in  a  dilute  solution  in 
equation  (4-23)  results  in 


T  d  19  9  R 

ip  (P  ,T)  +  kT  in  Cs  -  j  oo  r  ^  =  A 


(4-24) 


Differentiating  expression  (4-24)  with  respect  to  'r'  at 
constant  temperature  yields 


dip 


BP 


dP 
LR  dr 


LR 


+ 


kT  dCs 

C  dr 
s 


oo  rM2  =  0 


(4-25) 


j  A  ^  R 

where  2  -  q  '  s^nce  ^ 2  a  constant*  Substituting 


dr 


LR  . 


expression  (4-22)  for  dP  in  (4-25)  ,  and  rearranging  gives 

dr 


dC 


kT 


/llJf  L  dip  ,  2  , 

(M  -  p  )  03  rdr 

2  3PLR 


(4-26) 


From  equation  (4-12)  , 


dip 


By. 


BP 


LR 


BP 


LR 


Using  the  relations  y2 


3G_ 

BN. 


T,N1,N2 


and 


T,P,N. 


(4-27) 


V  = 


3G 

3P 
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[45]  reduces  equation  (4-27)  to 


t,n1,n2 


dip 


dp 


LR 


=  V, 


(4-28) 


where  is  the  partial  molar  volume  of  the  gas  in  the 
dilute  solution,  i.e. 


V  =  ^ 
2  3N. 


(4-29) 


T,P ,N. 


Substituting  expression  (4-29)  into  (4-26)  results  in 


dC  t _  9 

kT  =  (M2  -  P  V2)  oo  r  dr 


(4-30) 


Integrating  equation  (4-30)  from  r=Z  where  C  =  C  ,  to  a 

s  so 

point  ' r’  in  the  solution  gives 


C 

s 


C  exp 
so 


2 

oo 

2kT 


,  o  2  2  w  L- 

U  ~r  ) (p  V2 


(4-31) 


where  pL  and  V 2  are  assumed  to  be  constant.  The  concentra¬ 
tion  C  can  be  evaluated  using  Henry ' s  law  constant.  K 
so  2  H 

[46] ,  i.e. 


2 

The  symbol  G  denotes  the  Gibbs  potential,  which  is 
the  Legendre  transformation  of  the  internal  energy  replacing 
the  entropy  by  the  temperature  and  volume  by  the  pressure  as 
the  independent  variables,  i.e. 

G  =  U  -  TS  -  PV 
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C 

so 


P 

o 


VT) 


(4-32) 


once  Pq  is  specified.  Before  returning  to  the  primary 

system,  the  contribution  of  the  reference  system  can  be 

briefly  summarized.  Equation  (4-21)  represents  a  useful 

L  G 

substitution  for  the  unknown  quantity  ip  (P  ,T)  -  1-L02 

in  equation  (4-17)  for  the  primary  system.  The  reference 

concentration,  C  ,  expressed  by  equation  (4-31),  is  valid 

s 

for  the  case  of  negative  pressures  in  the  liquid-gas  solution 
outside  the  bubble. 

Returning  to  the  primary  system,  let  the  reference 
pressure  in  equation  (4-17)  be  equal  to  P  .  Then,  using 
equation  (4-21)  in  (4-17)  results  in 

pGx 

kT  Jin  (-=-^)  =  kT  Jin  (£-)  +  \  io2(r2  -  A2)M2  (4-33) 

o  s 

Substituting  equations  (4-31)  and  (4-32)  in  (4-33),  and 
rearranging  yields  an  expression  for  the  partial  pressure 
of  the  gas  inside  the  bubble,  i.e. 


PGX. 


C  K  exp 

n 


L  2w 

P  03  V2 

2kT 


(4-34) 


Using  equations  (4-16)  and  (4-34)  for  the  partial 
pressures  in  (4-8)  results  in  an  expression  for  the  total 
pressure,  PG,  inside  the  bubble  in  terms  of  the  properties 
of  the  liquid-gas  solution,  i.e. 


' 
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p  =  1lPv  +  ^2  C  KH 


where 


nn  =  exp 


—  (PL  -  P  )  -  C 
kT  v  v'  u 


and 


L  2- 

r  P  w  V 


n2  =  exP 


2kT 


*  (r2  -  *2) 


(4-35) 


(4-36) 


(4-37) 


Expression  (4-3)  for  'h'  requires  the  value  of 

PG  at  the  center  of  the  system,  PG.  Setting  ' r'  equal  to 

c 

zero  in  equation  (4-35)  gives 


pc  =  qpv  +  «2  C  KH  (4-38> 

where  and  are  defined  to  be  the  values  of  and  x\  , 
respectively,  at  r=0 ,  i.e. 


and 


5 


exp 


—  (PL  -  P  ) 
kT  v  c  v; 


-  C 


exp 


L  20 

p  co  Z  V, 

d 

2kT 


( 4-39 ) 


(4-40) 


Due  to  the  low  density  of  the  gas-vapor  mixture  inside  the 

Q 

bubble,  the  pressure  P  was  previously  assumed  to  be  constant 
and  equal  to  the  value  at  the  center.  It  should  be  noted 
that  this  is  equivalent  to  approximating  ti^  anc^  h2  ky  their 
values  at  the  center,  i.e. 
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n 


l 


(4-41) 


n2  =  (4-42) 

2  2 

where  r  <<  £  in  equations  (4-36)  and  (4-37). 3 

Equation  (4-40)  relates  £  to  the  angular  velocity 
and  the  length  of  the  system.  For  a  given  system  operating 
at  a  specific  speed,  both  03  and  £  are  constants.  Since  it 
is  preferable  to  express  as  a  function  of  the  thermo¬ 
dynamic  variables,  equation  (4-4)  will  be  used  to  eliminate 
1  L  2  2 

the  term  ^  p  co  £  from  equation  (4-40)  .  The  resultant 
equation  expresses  E,  as  a  function  of  the  pressure  in  the 
dilute  solution,  i.e. 


2 


exp 


(Pc  -  VV2 

kT 


(4-43) 


The  concentration  C  in  equation  (4-38)  specifically 
refers  to  that  in  the  liquid-gas  solution  at  the  center 
of  the  system.  In  practice,  the  solution  is  initially 
placed  in  the  tube  of  the  centrifugal  system  at  rest, 


3 

From  equation  (4-4)  ,  a  negative  pressure  in  the 
liquid-gas  solution  at  the  center  of  the  system  requires  that 


1  L  20  2 

J  P  03  £ 


> 


p 

o 


Due  to  practical  limitations  on  the  angular  velocity,  03,  in 
experimental  systems,  the  length,  £,  is  usually  several  orders 
of  magnitude  greater  than  the  equilibrium  size  of  the  bubble. 
Thus,  for  experimental  systems  such  as  that  of  Briggs  [47]/ 
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and  the  concentration  is  everywhere  the  same.  Then  the 

system  is  quickly  accelerated  to  the  angular  velocity 

required.  Since  the  time  scale  for  the  diffusion  of  a 

gas  in  a  liquid  is  much  greater  than  that  for  the  attain- 

4 

ment  of  mechanical  equilibrium,  the  concentration  of  the 
gas  remains  constant  throughout  the  solution  over  the  time 
span  of  the  experiment.  Thus,  C  is  considered  to  be  a  given 
quantity  in  expression  (4-38) . 

Substituting  equation  (4-38)  back  into  (4-3)  yields 
an  expression  for  'h'  in  terms  of  the  properties  of  the 
dilute  solution  for  a  given  system,  i.e. 


h 


2y 


'lPv  +  ?2  C  KH  -  p« 


(4-44) 


Briefly  summarizing,  for  an  equilibrium  gas-vapor 
bubble  subjected  to  a  centrifugal  field,  the  radius  of 
curvature,  'h',  at  the  center  of  the  system  is  the  basic 
parameter  characterizing  the  equilibrium  size  and  shape. 
For  the  special  case  of  a  dilute  solution  of  a  gas  in  a 
liquid,  the  introduction  of  a  reference  condition  valid 


^The  characteristic  time  scale 
process  is  of  the  order  of 

2 


TD  " 


D 


for  the  diffusion 


where  Dcis  the  diffusion  coefficient  of  the  dissolved  gas. 

For  a  solution  of  air  in  water  (Dcv  2.0  x  10  cm* 2  sec  1 * * *  [48]) 

with  l  equal  to  15  cm,  TD  is  of  the  order  of  107  seconds,  or 

about  3000  hours. 
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for  liquid  tensions  enables  'h'  to  be  expressed  in  terms 
of  the  properties  of  the  dilute  solution,  i.e.  T,  C,  P  , 
and  P^.  In  practice,  when  the  bubble  is  essentially 
spherical  '  h'  represents  the  critical  radius. 

In  the  next  chapter,  the  expression  developed  for 
' h'  will  be  used  to  calculate  the  critical  radius  of  a 
specific  system. 


CHAPTER  5 


THE  CRITICAL  RADIUS  OF  A  GAS-VAPOR  BUBBLE 
IN  A  DILUTE  SOLUTION  OF  AIR  IN  WATER 
AT  NEGATIVE  PRESSURES 

The  previous  chapter  treated  the  special  case  of  a 
gas-vapor  bubble  in  a  dilute  solution  of  a  single  gas  in 
a  liquid  solvent  at  negative  pressure.  In  practice  when 
e  <<  1,  the  radius  of  curvature  ' h'  at  the  center  represents 
the  critical  radius  of  the  bubble.  Equation  (4-44)  expresses 
' h'  as  a  function  of  the  properties  of  the  dilute  solution. 
Consider  now  the  application  of  the  expression  for  '  h'  to 
a  specific  system  containing  a  dilute  solution  of  air  in 
water  at  constant  temperature.'*'  Although  air  is  a  mixture 
of  primarily  two  gases,  it  is  commonly  treated  as  a  single 
gas  component  of  fixed  ratio,  i.e.  twenty-one  per  cent  oxygen 
and  seventy-nine  per  cent  nitrogen. 

Before  proceeding  to  evaluate  'h*  for  the  case  of 
air  in  water,  consider  the  range  of  validity  of  ' h' ,  the 
approximation  for  the  critical  radius.  As  the  bubble  tends 
to  a  spherical  shape,  the  ratio  of  the  maximum  radius,  'a', 
to  the  radius  of  curvature,  ' h',  at  the  center  approaches 

1Air  in  water  is  a  common  example  of  a  dilute 
solution  of  a  "gas"  in  a  liquid  for  which  an  abundance  of 
data  exists.  The  values  of  the  properties  of  the  system 
used  in  the  evaluation  of  ' h'  are  recorded  in  Appendix  C. 
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unity.  Equation  (2.2-14)  can  be  rearranged  to  give 


where 


a 

h 


1  +  f  (f)3 


L  2,  3 

P  03  h 

8y 


(5-1) 


(5-2) 


Figure  5-1  illustrates  the  deviation  of  the  ratio  —  from 
unity  as  a  function  of  the  parameter  ' f 1 .  Once  '  h'  is 
known,  the  value  of  'f'  follows  from  equation  (5-2).  For 
f  <<  1  in  Figure  5-1,  the  deviation  of  ^  from  unity  is 

3. 

approximately  equal  to  ' f  * .  As  —  approaches  unity,  *  f  * 
becomes  approximately  equal  to  the  shape  parameter  'e'.  Then 
for  f  <<  1,  i.e.  e  <<  1,  the  difference  between  'a'  and  'h' 
is  negligible,  and  the  value  of  'h'  represents  the  critical 
radius  of  the  bubble. 

Equation  (4-44)  expresses  'h'  as  a  function  of 

four  variables,  namely  the  temperature,  concentration, 

reference  pressure,  and  the  pressure  at  the  center  of  the 

system  in  the  dilute  solution.  The  variable  represents 

the  pressure  outside  the  bubble  at  the  center,  which  can  be 

positive  or  negative.  For  the  case  of  positive  pressures. 

Ward  et.  al.  [49]  expressed  the  critical  radius  as  a  function 

of  the  temperature,  concentration,  and  pressure  in  the  dilute 

solution.  At  zero  angular  velocity,  where  PL  is  equal  to 

c 

P  ,  expression  (4-44)  for  'h'  reduces  to  that  of  Ward  et  al . 


* 
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Deviation  of  the  Gas-Vapor  Bubble  From  Sphericity. 

r-  is  the  ratio  of  the  maximum  radius  to  the 
h 

radius  of  curvature  at  the  center;  f  is  the 
parameter  ^2^ 

8y 


Figure  5-1 
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When  the  system  is  rotating,  the  pressure  difference 

-  PQ  in  the  E, ^  term  in  expression  (4-43)  can  be  related 
to  the  angular  velocity  by  equation  (4-4).  Thus,  the 
additional  variable  Pq  in  the  expression  for  the  critical 
radius  is  a  result  of  the  pressure  variation  in  the  dilute 
solution  due  to  the  angular  rotation. 

In  presenting  the  results,  the  variation  of  'h'  with 

PL,  C,  and  T  will  be  examined,  since  these  variables  are 
c 

thermodynamic  properties  of  the  liquid  phase  at  the  center 
of  the  system.  However,  first  the  dependence  of  'h'  on  Pq 
must  be  determined.  The  reference  pressure  only  appears  in 
the  exponential  factor  E,^  equation  (4-44)  for  'h'.  From 
equation  (4-38),  the  term  £0CK  physically  represents  the 
increase  in  pressure  inside  the  bubble  due  to  the  presence 
of  the  solute  in  the  system. 

The  variation  of  'h'  with  Pq  is  presented  in  Figures 
5-2a,  5-2b,  and  5-2c  as  a  function  of  PL,  C,  and  T,  respect- 
ively.  Figure  5-2a  illustrates  the  effect  of  P^  on  the 
dependence  of  '  h'  on  Pq  for  constant  C  and  T.  The  change 
in  'h'  with  Pq  only  becomes  appreciable  as  P^  becomes 
positive,  and  near  the  upper  limit  considered  for  Pq. 
Physically,  the  effect  of  Pq  is  negligible  when  the  pressure 
outside  the  bubble  is  much  greater  in  magnitude  than  that 
inside . 

Next,  consider  the  effect  of  C  on  the  variation  of 

'  h'  with  P  for  constant  PL  and  T.  In  Figure  5-2b  the 
o  c 
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Figure  5-2a  Theoretical  Results  for  the  Dependence  of  the  Critical  Radius  on 

the  Reference  Pressure  at  Various  Negative  Pressures.  For  a 
solution  of  air  in  water  at  a  concentration  of  10“^  and  a  temperature 
of  25 °C 


concentration  of  dissolved  gas 
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Figure  5-2b  Theoretical  Results  for  the  Dependence  of  the  Critical  Radius  on 

the  Reference  Pressure  at  Various  Concentrations.  For  a  solution 
of  air  in  water  at  a  negative  pressure,  pf>  of  “10.13  MPa  and  a 
temperature  of  25°C 
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change  in  ' h'  with  P  increases  with  C.  This  is  evident 

o 

from  equation  (4-44)  where  the  relative  magnitude  of  the 
partial  pressure  of  the  gas,  i.e.  £_CK„,  in  comparison  to 
the  other  pressure  terms  in  the  denominator  determines  the 
significance  of  P  .  As  C  goes  to  zero,  the  effect  of  Pq 
disappears.  The  reason  the  curves  collapse  onto  the  pure 
vapor  limit  is  that  for  a  given  value  of  Pq  the  difference 
P^-Pq  becomes  increasingly  large  and  negative  as  Pq 
increases,  tending  to  cause  £  to  9°  to  zero.  The  maximum 
change  in  'h'  is  less  than  seven  per  cent  over  the  entire 
range  of  P  . 

Finally,  the  effect  of  T  on  the  change  in  ' h'  with 

P  at  constant  PL  and  C  is  presented  in  Figure  5-2c.  The 
o  c 

value  of  ' h'  for  a  given  temperature  is  almost  constant. 
Since  K  reaches  a  maximum  around  100°C,  the  partial  pressure 
of  the  gas  and  hence  the  effect  of  Pq  is  greatest  at  temper¬ 
atures  around  100°C.  Even  then,  the  change  in  ' h'  with  Pq 
is  only  perceptible  near  the  upper  limit  of  Pq  considered. 

Summarizing  the  results  above,  the  effect  of  P  on 

o 

1 h'  is  negligible  over  a  wide  range  of  values  of  Pq, 

especially  for  large  negative  pressures.  Since  'h'  is 

independent  of  P  ,  the  critical  radius  can  be  considered  a 

function  of  PL,  C,  and  T.  The  value  of  P  used  for  computa- 

c  o 

tional  purposes  will  be  the  equilibrium  pressure  of  a  gas- 
vapor  mixture  over  a  flat  interface  of  a  dilute  solution  at 
temperature  T  and  concentration  C.  Then,  PQ  is  expressed  by 
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the  relation  [50] 


where 


=  n  P  +  C  K 
o  v  H 


0 


.  v 


exp  [kT  (Po  ' 


p  )  -  C] 
v 


(5-3) 


(5-4) 


Experimentally,  this  value  of  Pq  corresponds  to  the  pressure 

at  which  the  liquid-gas  solution  is  prepared.  Since  P^  and 

Ktt  depend  on  T,  P  is  a  function  of  T  and  C.  However,  since 
H  o 

' h'  is  independent  of  Pq  over  the  range  of  Pq  considered, 
the  effect  on  ' h'  of  the  change  in  Pq  with  T  and  C  is 
negligible . 

It  should  be  noted  that  the  choice  of  Pq  given  by 

equation  (5-3)  is  equal  to  the  pressure  inside  the  bubble  at 

zero  angular  velocity.  Since  the  pressures  inside  and  outside 

the  bubble  are  then  equal,  the  equilibrium  interface  for 

zero  angular  velocity  has  an  infinite  radius  of  curvature, 

i.e.  is  a  flat  surface.  Thus,  as  PL  approaches  the  value  of 

c 

Pq  given  by  equation  (5-3),  the  critical  radius  tends  to 
infinity,  since  the  denominator  in  equation  (4-44)  approaches 


zero . 


Having  established  the  value  of  Pq  to  be  used,  the 

dependence  of  'h'  on  P^,  C,  and  T  can  now  be  investigated. 

The  variation  of  ' h '  with  PL  for  several  values  of  C  and  T 

c 

is  presented  in  Figures  5-3a  to  5-3d.  The  value  of  ' h* 

decreases  as  the  magnitude  of  the  negative  pressure  PL 

o 


Critical  Radius,  h  (pm) 
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Negative  Pressure,  (MPa) 


Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Negative  Pressure  at  Various 
Concentrations.  For  a  solution  of  air  in  water 
at  a  temperature  of  25°C 


Figure  5-3a 


Critical  Radius,  h  (yin) 
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Negative  Pressure,  PL  (MPa) 

c 


Figure  5-3b  Theoretical  Results  for  the  Dependence  of  the 

Critical  Radius  on  Negative  Pressure  at  Various 
Concentrations.  For  a  solution  of  air  in  water 
at  a  temperature  of  100°C 
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Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Negative  Pressure  at  Various 
Concentrations.  For  a  solution  of  air  in  water 
at  a  temperature  of  200 °C 


Figure  5-3c 
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Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Negative  Pressure  at  Various 
Concentrations.  For  a  solution  of  air  in  water 
at  a  temperature  of  300 °C 


Figure  5-3d 
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increases.  In  equation  (4-4),  as  P  becomes  large  and 

c 

negative,  the  pressure  difference  across  the  interface 

steadily  increases,  causing  the  critical  radius  to  decrease. 

The  rate  of  change  of  ' h'  with  PL  increases  as  PL  becomes 

c  c 

positive.  The  value  of  ' h'  tends  to  infinity  as  the  pressure 

in  the  liquid  tends  to  the  value  P  . 

o 

Figures  5-3a  to  5-3d  also  illustrate  the  effect  of 

the  solute  concentration  on  the  critical  radius.  At  constant 

PL  and  T,  the  value  of  '  h'  decreases  as  C  increases.  For 
c 

PL  large  and  negative,  only  high  concentration  levels 
c 

appreciably  reduce  the  value  of  'h'.  However,  the  effect  of 
C  on  1  h'  increases  as  P^  becomes  positive. 

The  effect  of  increasing  C  is  to  increase  the  partial 
pressure  of  the  gas  inside  the  bubble,  thus  increasing  the 
pressure  difference  across  the  bubble  interface,  and  decreas¬ 
ing  the  critical  radius.  As  P^  becomes  smaller  in  magnitude, 
the  significance  of  the  partial  pressure  of  the  gas  increases, 
and  the  curves  for  each  concentration  tend  to  diverge. 

Comparing  the  spread  of  the  curves  in  each  of  Figures 
5-3a  to  5-3d  indicates  that  the  dependence  of  the  value  of 
' h'  on  C  decreases  with  T.  Since  K„  decreases  for  temperatures 

tl 

above  100°C,  the  partial  pressure  of  the  gas  also  decreases 

with  T.  This  explains  the  reduced  effect  of  C  on  'h'  with  T. 

The  results  are  summarized  in  Figure  5-4  which 

illustrates  the  change  in  ' h'  with  C  at  constant  PL  and  T. 

c 

-3 

An  upper  bound  of  10  has  been  used  for  C  in  order  to  better 
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illustrate  the  effect  of  C  on  ' h ' ,  although  this  value  is 

bordering  on  the  upper  limit  for  a  dilute  solution  [51] . 

For  C  close  to  zero,  the  curves  for  negative  values  of  PL 

c 

are  almost  horizontal.  The  exact  magnitudes  of  the  change  in 
' h'  with  C  for  these  values  are  presented  in  Table  1. 

The  last  effect  to  be  considered  is  that  of  tempera¬ 
ture  on  'h'  at  constant  PL  and  C.  In  equation  (4-44)  the 

o 

three  properties  which  change  significantly  with  T  are  the 
surface  tension,  vapor  pressure,  and  Henry's  law  constant. 
Table  2  shows  that  y  decreases  with  T,  while  P  increases. 

The  value  of  K^.  initially  increases  up  to  a  maximum  between 

ri 

75°C  and  100°C,  and  thereafter  steadily  decreases.  The 

solubility  has  a  minimum  value  at  the  maximum  value  of  K  . 

H 

The  temperature  variation  of  ' h'  is  illustrated  in 

Figures  5-5a  and  5-5b  for  the  pure  vapor  case  and  a  concen- 

-3 

tration  of  10  ,  respectively.  The  pure  vapor  case  repre¬ 

sents  an  upper  limit  for  ' h'  at  constant  P^  and  T.  The 
temperature  variation  of  'h'  for  the  pure  vapor  case  depends 
only  on  the  change  in  y  and  P^  with  T.  From  Figure  5-5a, 

the  value  of  'h'  decreases  with  T  for  constant  PL.  The  rate 

c 

of  change  of  'h'  with  T  increases  with  temperature,  and  as 
P^  becomes  positive,  since  the  relative  magnitude  of  P^ 
compared  to  P^  becomes  larger. 

The  same  general  decrease  in  'h'  with  temperature 

-3 

is  shown  in  Figure  5-5b  for  a  concentration  of  10  .  For  a 

dilute  solution,  the  temperature  variation  of  'h'  depends 
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TABLE  1 

THE  DEPENDENCE  OF  THE  CRITICAL  RADIUS  ON  THE  AIR  CONTENT 
FOR  WATER  WITH  VARIOUS  NEGATIVE  PRESSURES 
AT  A  TEMPERATURE  OF  25 °C 


Dissolved 

Gas 

Content 

C 

(mole  ratio) 

-10.13 

Critical 

Negative 

-5.07 

Radius,  h 

Pressure ; 

-2.53 

(ym) 

PL  (MPa) 
c 

-0.10 

o 

• 

o 

1  x 

10-3 

0 . 00901 

0.01272 

0.01560 

0.01891 

0.02152 

1  X 

1 

o 

i — i 

0.01336 

0.02503 

0.04440 

0 . 08288 

0.3889 

1  X 

10"5 

0.01411 

0 .02802 

0.05521 

0.01322 

1.873 

1  X 

VD 

1 

O 

i — 1 

0.01419 

0.02836 

0.05660 

0.1406 

13.64 

1  X 

O 

1 

'-J 

0 . 01420 

0.02839 

0.05674 

0.1415 

36.75 

0 

0.01420 

0.02840 

0.05676 

0.1416 

45.28 
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Figure  5-5a  Theoretical  Results  for  the  Dependence  of  the 

Critical  Radius  on  Temperature  for  a  Pure  Vapor 
Bubble  in  Water  at  Various  Negative  Pressures 


Critical  Radius,  h  (ym) 


108 


Theoretical  Results  for  the  Dependence  of  the 
Critical  Radius  on  Temperature  at  Various  Negative 
Pressures.  For  a  solution  of  air  in  water  at  a 
concentration  of  10”^ 


Figure  5-5b 
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upon  the  change  in  K„  with  T,  in  addition  to  the  change  in 

H 

y  and  P  .  The  relative  variation  of  the  magnitudes  of  the 

partial  pressures  of  the  vapor  and  gas  inside  the  bubble 

due  to  the  temperature  variation  of  P^  and  introduces- 

an  inflection  point  in  the  curves  of  constant  pressure  in 

Figure  5-5b.  The  inflection  point  occurs  around  100°C,  at 

which  temperature  the  value  of  K„  undergoes  a  maximum.  As 

H. 

PL  becomes  positive,  the  temperature  variation  of  the  partial 
o 

pressures  becomes  more  significant.  In  particular,  the 

decrease  in  K  tends  to  offset  the  temperature  variation  of 

'h'  with  y  and  Py  between  100°C  and  200°C.  For  large 

temperatures,  the  effect  of  P^  dominates,  reducing  the  value 

of  'h'  sharply.  The  dependence  of  ' h'  on  the  relative 

variation  of  the  partial  pressures  becomes  less  and  less 

as  the  magnitude  of  the  negative  pressure  increases,  since 

only  the  pressure  difference  across  the  interface  is  important. 

The  value  of  the  critical  radius  at  the  critical 

temperature  represents  a  lower  limit  for  the  temperature 

variation  of  'h'.  At  the  critical  point  for  the  liquid-gas 
2 

solution,  the  liquid  and  vapor  phases  become  identical. 

Then,  the  interface  between  them  ceases  to  exist,  and  the 
surface  tension  becomes  equal  to  zero.  Thus,  the  value  of 

o 

The  values  of  T  and  P  for  the  dilute  solution  of 

air  in  water  are  not  equal  to  tfie  values  for  pure  water,  i.e. 

T  =  374°C  and  P  =  22.12  MPa  [52].  In  practice,  a  "mixing 
rule"  is  used  torelate  the  value  of  T  and  Pc  for  the  solution 
to  the  values  of  T  and  P  for  the  pure  solvent  and  solute 

r  r-  n  C  C 
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the  critical  radius  given  by  equation  (4-44)  is  equal  to 

zero  at  the  critical  point. 

The  value  of  ' h'  is  inversely  proportional  to  the 

pressure  difference  across  the  bubble  interface.  The 

expressions  in  the  denominator  of  equation  (4-44)  for  the 

partial  pressures  of  the  gas-vapor  mixture  inside  the  bubble 

include  the  factors  E, and  E,^,  given  by  expressions  (4-39) 

and  (4-43),  respectively.  These  factors  are  functions  of 

PL,  C,  T,  and  P  .  Since  the  choice  of  P  given  by  equation 
c  o  o 

(5-3)  is  a  function  of  C  and  T,  it  is  sufficient  to  consider 

the  variation  with  PL,  C,  and  T. 

c 

The  exponential  terms  E,^  and  E,  represent  reductions 
in  the  partial  pressures  of  the  vapor  and  gas  inside  the 
bubble  due  to  the  curvature  of  the  interface.  From  equation 
(5-3) ,  for  a  gas-vapor  mixture  in  equilibrium  with  a  dilute 
solution  across  a  flat  interface,  the  partial  pressures  of 
the  vapor  and  gas  are  given  respectively  by 

p°  =  n  pv  (5-5) 

and 

P2  =  C  Kh  (5-6) 

3 

For  n  approximately  equal  to  one,  i.e.  the  partial  pressure 
of  the  vapor  approximately  equal  to  the  vapor  pressure,  the 


of 


For  a  dilute  solution 
25°C  and  a  concentration  of 


of  air  in 


10 


“4 


the 


water  at  a  temperature 
value  of  n  is  1.0053. 


Ill 


percentage  reductions  in  the  partial  pressures  of  the  vapor 

and  gas  due  to  the  curvature  of  the  interface  become  (1  -  E,^) 

and  (1  -  E,  ) ,  respectively. 

Table  3a  illustrates  the  variation  of  (1  -  £^)  anc^ 

(1  -  E,^)  with  for  constant  C  and  T.  Each  value  increases 

steadily  as  the  magnitude  of  the  negative  pressure  increases. 

The  deviation  of  and  E,^  from  unity  for  large  negative 

pressures  significantly  lowers  the  pressure  inside  the  bubble 

in  comparison  to  that  for  the  case  of  a  flat  interface. 

The  variation  of  (1  -  and  (1  -  E,^)  with  C  at 

constant  PL  and  T  is  presented  in  Table  3b.  The  value  of 
o 

(1  -  E,^)  is  almost  constant,  while  that  of  (1  -  E,  )  decreases 

as  c  increases.  In  expression  (4-39)  for  £  ,  the  effect  of 

C  is  important  only  when  P^  approaches  P  .  In  expression 

(4-43)  for  £  ,  the  only  effect  of  C  is  to  change  P  .  This 

effect  is  negligible  when  the  magnitude  of  PL  is  much  greater 

c 

than  P  . 
o 

Finally,  Table  3c  considers  the  temperature  variation 

of  (1  -  £,)  and  (1  -  E,  )  for  constant  PL  and  C.  The  values 
X  ^  c 

of  (1  -  £^)  and  (1  -  E,^)  initially  decrease  to  a  minimum, 

and  then  increase  with  T ,  as  a  result  of  the  temperature 

variation  of  vL,  P  ,  and  P  ,  and  the  presence  of  T  in  the 

denominator  of  the  exponents  of  both  E,  and  £  .  The  exact 

dependence  of  the  percentage  reductions  in  the  partial 

pressures  represented  by  (1  -  £^)  and  (1  -  Ea^)  on  both  C  and 

T  is  sensitive  to  the  value  of  PL.  However,  in  general  the 

c 


112 


TABLE  3a 

THE  PERCENTAGE  REDUCTIONS  IN  THE  PARTIAL  PRESSURES 
FOR  VARIOUS  NEGATIVE  PRESSURES 
(For  a  solution  of  air  in  water  with  a  concentration  of  10 

and  a  temperature  of  100°C) 


Reduction  in  Partial  Pressures 
Vapor  Gas 

1  -  C-l  (%)  1  *  Z2  (%) 


o 

• 

o 

0.072 

1.239 

o 

• 

rH 

i 

0.675 

2.284 

-2.5 

1.573 

3.831 

-5.0 

3.051 

6.354 

-7.5 

4.507 

8.812 

o 

• 

o 

1 — 1 

1 

5.942 

11.20 

-25.0 

14.11 

24.30 

-50.0 

26.17 

41.98 

-75.0 

36.54 

55.52 

-100.0 

45.45 

65.91 

Negative 

Pressure 

PL  (MPa) 
c 


113 


TABLE  3b 

THE  PERCENTAGE  REDUCTIONS  IN  THE  PARTIAL  PRESSURES 

FOR  VARIOUS  CONCENTRATIONS 

(For  a  solution  of  air  in  water  with  a  negative  pressure,  PL, 

c 

of  -10  MPa  and  a  temperature  of  100°C) 


Concentration  of 
Dissolved  Gas 

C  (mole  ratio) 

Reduction 

Vapor 

1  -  (%) 

in  Partial  Pressure 

Gas 

1-C2  (%) 

1 

X 

10-3 

6.032 

19.87 

5 

X 

io"4 

5.985 

15.17 

1 

X 

io"4 

5.942 

11.21 

5 

X 

10-5 

5.938 

10.70 

1 

X 

io"5 

5.933 

10.29 

5 

X 

io-6 

5.933 

10.24 

1 

X 

io-6 

5.933 

10.20 
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TABLE  3c 

THE  PERCENTAGE  REDUCTIONS  IN  THE  PARTIAL  PRESSURES 

FOR  VARIOUS  TEMPERATURES 

(For  a  solution  of  air  in  water  with  a  negative  pressure,  PL, 

c 

-4 

of  -10  MPa  and  a  concentration  of  10  ) 


Temperature 

T  (°C) 

Reduction 

Vapor 

1  -  c1  <%) 

in  Partial  Pressures 

Gas 

1  -  s2  (%) 

50 

6.573 

12.61 

100 

5.942 

11.21 

150 

5.688 

10.10 

200 

5.942 

9.702 

250 

6 . 995 

10.32 

300 

9.404 

12.21 
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variation  of  (1  -  £  )  and  (1  -  with  C  and  T  is  small 

compared  to  the  variation  with  P^. 

Having  obtained  actual  values  for  the  critical 

radius  using  expression  (4-44)  for  'h',  the  validity  of 

this  approximation  can  now  be  determined.  Recall  that  the 

value  of  ' h'  for  the  pure  vapor  case  represents  an  upper 

limit  at  a  given  PL  and  T.  Then,  from  Figure  5-5a,  it  is 

c 

evident  that  for  PL  less  than  -5.0  MPa,  'h'  is  always  less 

c 

than  0.1  ym.  Using  this  value  of  'h'  in  equation  (5-2)  gives 

-12 

a  value  of  2  x  10  for  'f'  at  a  temperature  of  25°C  and 

an  angular  velocity  of  1000  rad/sec.  From  Figure  5-1,  the 

deviation  of  ^  from  unity  is  approximately  equal  to  the  value 

-12 

of  'f' ,  2  x  10  .  Since  this  deviation  can  be  considered 

negligible,  the  bubble  is  essentially  spherical  of  radius  'h' . 

In  summary,  using  expression  (4-44)  for  ‘h1,  the 

critical  radius  was  evaluated  for  a  dilute  solution  of  air 

in  water.  Since  'h'  is  insensitive  to  P  ,  the  critical 

o 

radius  is  essentially  a  function  of  PL,  C,  and  T.  The  value 

c 

of  ' h*  decreases  as  the  magnitude  of  the  negative  pressure 
in  the  solution  outside  the  bubble  increases.  Increasing 
the  solute  concentration  decreases  the  value  of  1 h'.  Finally, 
' h'  decreases  significantly  at  negative  pressures  as  the 
temperature  increases.  In  the  region  of  negative  pressures 
the  bubble  is  essentially  spherical. 
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CHAPTER  6 


SUMMARY  AND  CONCLUSIONS 

In  the  beginning  of  this  investigation  the  concept 
of  a  critical  size  gas-vapor  bubble  was  introduced  to  explain 
the  phase  transition  of  a  liquid-gas  solution  to  a  gas-vapor 
mixture.  Homogeneous  nucleation  theory  postulates  the 
creation  of  minute  gas-vapor  nuclei  by  fluctuations  in  the 
density  of  the  liquid-gas  solution.  The  reversible  work 
required  to  form  a  critical  size  nucleus  represents  a  type  of 
"potential  barrier"  to  the  transition  of  the  liquid-gas 
solution  to  a  gas-vapor  mixture  [54] .  If  a  bubble  smaller 
than  the  critical  size  is  formed,  it  spontaneously  decays. 

If  a  bubble  larger  than  the  critical  size  is  formed,  it 
spontaneously  grows,  becoming  a  center  for  the  formation  of 
the  new  phase.  Thus,  the  idea  of  a  critical  size  nucleus  is 
central  to  homogeneous  nucleation  theory. 

For  the  case  of  positive  pressure  in  the  solution, 
the  bubble  is  spherical,  and  the  critical  radius  has  pre¬ 
viously  been  determined.  However,  a  solution  can  be  in  a 
metastable  state  at  negative  pressures,  i.e.  the  solution  can 
be  in  tension.  The  purpose  of  this  investigation  is  to 
determine  the  critical  size  of  a  gas-vapor  bubble  for  the 
case  of  a  solution  at  negative  pressure. 
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The  thermodynamic  system  considered  corresponds  to 
that  used  in  practice  to  measure  the  tensile  strengths  of 
liquids  by  means  of  a  centrifugal  field.  The  gas-vapor 
bubble  is  imagined  to  exist  at  the  center  of  rotation  of  a 
narrow  tube  filled  with  a  liquid-gas  solution.  If  the  tube 
is  rotated  with  sufficient  angular  velocity,  the  liquid 
outside  the  bubble  at  the  center  of  the  system  can  be  placed 
in  a  state  of  tension. 

The  critical  size  and  shape  of  the  equilibrium  bubble 
follow  from  the  conditions  of  thermodynamic  equilibrium.  A 
plane  section  through  the  axis  of  rotation  of  the  bubble 
can  be  approximated  for  negative  pressure  by  the  ellipse 
given  by 

2  2 

- ~  2  2  +  H  =  1  (6-D 

(1  +  e)  a  a 

where  'a'  is  the  maximum  radius  of  the  bubble,  and  '  e'  is 
a  shape  parameter  defined  by  equation  (2.2-15).  The 
magnitude  of  ' e'  determines  the  deviation  of  the  bubble 
from  sphericity.  For  all  practical  purposes,  at  negative 
pressures  the  magnitude  of  ' e'  is  negligible,  which  implies 
that  the  bubble  is  essentially  spherical. 

The  equilibrium  state  of  the  system  is  an  unstable 
equilibrium,  since  the  free  energy  of  the  bubble  is  a 
maximum.  Thus,  the  equilibrium  size  of  the  bubble  represents 
a  critical  value.  A  smaller  bubble  will  spontaneously 


. 
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collapse,  while  a  larger  bubble  will  spontaneously  expand. 

A  thermodynamic  fluctuation  in  the  density  of  the  solution 
which  produces  a  bubble  larger  than  the  critical  size  will 
initiate  the  formation  of  a  gas-vapor  filled  cavity.  The 
reversible  work  required  to  form  a  critical  size  gas-vapor 
bubble  is  given  by 

WR  =  'I’l-Y-J*2  (6-2) 

where  for  negative  pressures  the  magnitude  of  'e'  is 
negligible . 

The  critical  radius  is  closely  approximated  by  the 
radius  of  curvature,  '  h* ,  at  the  center  of  the  bubble,  since 
the  value  of  'e'  is  very  small.  Thus,  the  critical  radius 
can  be  represented  by  the  expression  for  ' h'  obtained  from 
the  equilibrium  conditions.  For  the  special  case  of  a  dilute 
solution  of  a  gas  in  a  liquid,  the  critical  radius  can  be 
expressed  in  terms  of  the  properties  of  the  dilute  solution. 
Introducing  an  appropriate  reference  condition  valid  for 
negative  pressures  in  the  dilute  solution  outside  the  bubble 
results  in 

h  =  - — — - r  (6-3) 

qpv +  52ckh  -  Pc 

where  and  are  exponential  terms  which  deviate  from 
unity,  and  are  expressed  by  equations  (4-39)  and  (4-43)  , 
respectively . 
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Using  equation  (6-3)  for  'h',  the  critical  radius 
was  evaluated  for  the  specific  case  of  a  dilute  solution  of 
air  in  water.  The  critical  radius  was  essentially  a  function 
of  the  temperature,  concentration,  and  pressure  of  the  liquid 
phase  at  the  center  of  the  system,  being  independent  of  the 
reference  pressure.  As  the  magnitude  of  the  negative 
pressure  increased,  the  value  of  the  critical  radius  decreased 
The  effect  of  the  solute  was  to  reduce  the  bubble  size, 
although  a  significant  reduction  was  only  effected  by  concen¬ 
trations  near  the  upper  limit  for  a  dilute  solution.  The 
dependence  of  the  critical  radius  on  the  temperature  followed 
from  the  variation  of  the  surface  tension,  vapor  pressure, 
and  Henry's  law  constant  with  temperature. 

The  dependence  of  the  expression  for  the  critical 
radius  given  by  (6-3)  on  the  solution  properties  outside  the 
bubble  has  several  implications  for  the  tensile  strength  of 
a  dilute  solution  in  a  centrifugal  system.  These  implica¬ 
tions  assume  that  the  liquid  fracture  occurs  in  the  body  of 
the  fluid,  i.e.  not  at  a  surface,  and  that  homogeneous  nuclea- 
tion  is  the  mechanism  responsible  for  the  formation  of  the  gas 
vapor  cavity  precipitating  liquid  fracture.  As  the  magnitude 
of  the  negative  pressure  in  the  solution  at  the  center  in¬ 
creases,  the  value  of  the  critical  radius  steadily  decreases. 
The  reversible  work  required  to  form  the  critical  nucleus 
decreases  as  the  critical  radius  becomes  smaller.  However, 
as  the  reversible  work  decreases,  the  probability  of  the 
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formation  of  a  critical  size  gas-vapor  nucleus  increases, 
since  the  probability  is  proportional  to  the  negative 
exponential  of  the  reversible  work.  At  some  point,  the 
formation  of  a  nucleus  greater  than  critical  size  becomes 
reasonably  likely  over  a  short  period  of  time.  When  this 
negative  pressure  is  reached,  the  liquid  at  the  center  would 
be  expected  to  rupture.  This  negative  pressure  represents 
the  tensile  limit  of  the  liquid  solution. 

The  fact  that  increasing  the  concentration  of  the 
gas  does  not  significantly  reduce  the  critical  radius  except 
at  large  concentration  levels  implies  that  for  a  dilute 
solution  the  tensile  strength  is  not  substantially  reduced 
by  the  presence  of  dissolved  gas.  Since  there  is  no  general 
concensus  as  to  the  effect  of  dissolved  gas  on  the  experi¬ 
mentally  determined  tensile  strengths  of  liquids,  the 
validity  of  this  implication  cannot  at  present  be  confirmed. 

The  temperature  dependence  of  the  critical  radius 
suggests  that  the  tensile  strength  of  a  liquid-gas  solution 
is  significantly  reduced  as  the  temperature  of  the  system 
increases.  Although  there  has  been  both  experimental  and 
theoretical  work  which  seems  to  bear  this  out  [55] ,  there  is 
a  notable  lack  of  data  from  a  centrifugal  system  at  high 
temperatures . 

This  investigation  considered  the  equilibrium  size 
and  shape  of  a- gas-vapor  bubble  in  a  liquid-gas  solution 
subject  to  a  centrifugal  field.  The  radial  pressure 
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distribution  in  the  solution  due  to  the  centrifugal  field 
was  used  to  generate  tensions  in  the  solution  outside  the 
bubble.  For  the  small,  equilibrium  gas-vapor  bubbles 
corresponding  to  negative  pressures  in  the  liquid-gas 
solution,  the  bubble  shape  is  approximately  spherical.  For 
such  a  bubble  in  a  dilute,  liquid-gas  solution,  the  critical 
radius  is  a  function  of  the  properties  of  the  liquid  phase 
outside  the  center  of  the  bubble.  Since  the  equilibrium 
size  of  the  gas-vapor  bubble  is  very  small,  any  direct 
observation  of  the  equilibrium  size  and  shape  would  be 
impossible.  However,  the  expression  for  the  critical  radius 
can  be  used  to  obtain  theoretical  predictions  for  the  tensile 
strengths  of  dilute,  liquid-gas  solutions  in  a  centrifugal 
system  on  the  basis  of  homogeneous  nucleation  theory.  The 
comparison  of  these  predictions  with  experimentally  deter¬ 
mined  values  would  provide  a  means  of  indirectly  verifying 
the  results  obtained  in  this  investigation. 


. 


v 


BIBLIOGRAPHY 


122 


BIBLIOGRAPHY 


Knapp,  R.T.,  Daily,  J.W.,  and  Hammitt,  F.G., 

"Cavitation",  McGraw-Hill,  New  York,  1970, 

p .  1 . 

Plesset,  M.S.,  "The  Tensile  Strength  of  Liquids", 
Cavitation  State  of  Knowledge,  ASME,  1969, 
p .  15 . 

Fischer,  J.C.,  "The  Fracture  of  Liquids",  J.  Appl. 
Phys.,  19,  1948,  p.  1064. 

Holl,  J.W.,  "Nuclei  and  Cavitation",  Trans.  ASME, 

J.  Basic  Eng.,  9_2,  Dec.  1970  ,  p.  681. 

Knapp,  R.T.  et  al. ,  op.  cit . ,  p.  62. 

Frenkel,  J.,  "Kinetic  Theory  of  Liquids",  Dover, 
New  York,  1955,  pp .  366-381. 

Landau,  L.D.,  and  Lifshitz,  E.M.,  "Statistical 
Physics",  Pergamon,  Oxford,  2nd  ed . ,  1968, 
pp.  471-474. 

Callen,  H.B.,  "Thermodynamics",  Wiley,  New  York, 
1960,  pp.  267-272. 

Landau,  L.D.,  and  Lifshitz,  E.M.,  op.  cit. ,  p.  349. 


Ward,  C.A.,  Balakrishnan ,  A.,  and  Hooper,  F.C., 

"On  the  Thermodynamics  of  Nucleation  in  Weak 
Gas-Liquid  Solutions",  Trans.  ASME,  J.  Basic 
Eng.,  92_,  Dec.  1970,  pp .  695-701. 

Blake,  F.G.,  Jr.,  "The  Tensile  Strength  of  Liquids: 

A  Review  of  the  Literature",  Tech.  Memo.  No.  9, 
Acoustics  Res.  Lab.,  Harvard  Univ. ,  June  11, 
1949  . 

Galloway,  W. ,  "An  Experimental  Study  of  Acoustically 
Induced  Cavitation  in  Liquids",  J.  Acoust.  Soc . 
Am.,  26,  1954,  pp.  849-857. 


Briggs,  L. ,  "Limiting  Negative  Pressure  of  Water", 
J.  Appl.  Phys.,  21,  1950,  pp .  721-722. 


124 


14.  Reynolds,  0.  ,  "On  the  Internal  Cohesion  of  Liquids 

and  the  Suspension  of  a  Column  of  Mercury  to 
a  Height  of  More  Than  Double  That  of  the 
Barometer",  Mem.  Manch.  Lit.  Phil.  Soc.  ,  1_, 
Ser.  3,  1882,  pp.  1-19. 

15.  Temperley,  H.N.V.,  and  Chambers,  LL.G. ,  "The 

Behaviour  of  Water  Under  Hydrostatic  Tension: 
I",  Proc.  Phys .  Soc.  (London),  5j3_,  1946  , 
pp.  420-436. 

16.  Apfel,  R.E.,  and  Smith,  M.P.,  Jr.,  "The  Tensile 

Strength  of  Di-ethyl  Ether  Using  Briggs's 
Method",  J.  Appl .  Phys.,  4_8_,  1977  ,  pp .  2077- 
2078  . 

17.  Donoghue,  J.J.,  Vollrath,  R.E.,  and  Gerjuoy,  E., 

"The  Tensile  Strength  of  Benzene",  J.  Chem. 
Phys.,  19,  1951,  pp .  55-61. 


18.  Hammitt,  F.G.,  "Effects  of  Gas  Content  Upon  Cavitation 

Inception,  Performance,  and  Damage",  Report  No. 
UMICH  01357-19-T,  Cavitation  and  Multiphase  Flow 
Lab.,  Univ.  of  Michigan,  1971,  p.  20. 

19.  Iyengar,  K.S.,  and  Richardson,  E.G.,  "The  Role  of 

Cavitation  Nuclei",  Fluids  Report  No.  57, 

Mech.  Engr.  Res.  Lab.,  Fluid  Mech.  Div., 

August,  1957. 

20.  Kuper,  C.G.,  and  Trevena,  D.H.,  "The  Effect  of 

Dissolved  Gases  on  the  Tensile  Strength  of 
Liquids",  Proc.  Phys.  Soc.  (London),  6_5,  1952, 
pp .  46-54 . 

21.  Harvey,  E.N.  et  a]^.  ,  "Bubble  Formation  in  Animals  II", 

J.  Cell.  Comp.  Physiol.,  _24_,  1944  ,  pp.  23-27. 

22.  Rayleigh,  Lord,  "The  Equilibrium  of  Revolving  Liquid 

Under  Capillary  Force",  Phil.  Mag.,  2_8,  1914, 
pp .  161-170. 

23.  Rosenthal,  D.K.,  "The  Shape  and  Stability  of  a  Bubble 

at  the  Axis  of  a  Rotating  Liquid",  J.  Fluid 
Mech.,  12_,  1962,  pp.  358-363. 

24.  Princen,  H.M.,  Zia,  I.Y.Z.,  and  Mason,  S.G., 

"Measurement  of  Interfacial  Tension  from  the 
Shape  of  a  Rotating  Drop",  J.  Coll.  Int.  Sci., 

23,  1967,  pp.  99-107. 


125 


25. 

Callen,  H.B.,  op.  cit.,  p.  105. 

26  . 

Gibbs,  J.W.,  "The  Scientific  Papers 
Gibbs",  Dover,  New  York,  1961, 

of  J.  Willard 
pp .  219-231. 

27. 

Ibid . ,  p .  227 . 

28  . 

Ibid . ,  p .  229 . 

29. 

Ibid.,  pp.  276-285. 

30  . 

Defay,  R. ,  Prigogine,  I.,  and  Bellemans,  A.,  "Surface 

Tension  and  Adsorption",  Wiley,  New  York,  1966, 
pp.  87-88. 


31. 

Kraus ,  H . , 

"Thin 

Elastic 

Shells",  Wiley,  New  York, 

1967  , 

p .  20 

• 

32. 

Rayleigh , 

Lord , 

op.  cit. 

,  p .  163 . 

33. 

Frenkel,  J 

.  ,  op. 

cit. ,  p 

.  374  . 

34. 

Daugherty , 

R  •  L  •  / 

and  Franzini,  J.B.,  "Fluid 

Mechanics  With  Engineering  Applications", 
McGraw-Hill,  New  York,  7th  ed.,  1977,  p.  540. 


35.  Spiegel,  M.R. ,  "Mathematical  Handbook  of  Formulas 

and  Tables",  Schaum's  Outline  Series,  McGraw- 
Hill,  New  York,  1968,  p.  110. 

36.  Callen,  H.B.,  op.  cit. ,  pp .  131-145. 


37.  Glansdorf,  P.,  and  Prigogine,  I.,  "Thermodynamic 

Theory  of  Structure,  Stability,  and  Fluctuations", 
Wiley,  London,  1971,  pp .  14-16. 


u> 

00 

• 

Ward . ,  C . A.  et 

al . ,  op .  cit . , 

p.  697. 

39. 

Ibid . ,  p .  696 . 

40. 

Landau ,  L . D . , 

and  Lifshitz,  E 

ol 

• 

£ 

• 

cit . ,  p . 

275. 

41. 

Prigogine ,  I . , 

and  Defay ,  R. , 

"Chemical  Thermodynamic 

Longmans 

Green,  London, 

1954,  pp 

.  126-127 

• 

42. 

Landau ,  L . D . , 

and  Lifshitz,  E 

.M. ,  op. 

cit . ,  p . 

277. 

Andrews,  F.C.,  "Thermodynamics:  Principles  and 

Applications",  Wiley,  New  York,  1971,  p.  162. 


43. 


126 


44.  Ward,  C.A.  et  al . ,  op.  cit. ,  p.  697. 

45.  Callen,  H.B.,  op.  cit. ,  pp .  99-100. 

46.  Prigogine,  I.,  and  Defay,  R. ,  op.  cit. ,  p.  334. 

47.  Briggs,  L. ,  op.  cit. ,  p.  721. 

48.  Bolz,  R.E.,  and  Tuve ,  G.L.,  eds . ,  "CRC  Handbook 

of  Tables  for  Applied  Engineering  Science", 

The  Chemical  Rubber  Co.,  Cleveland,  1970,  p.  453. 

49.  Ward,  C.A.  et  al . ,  op .  cit. ,  p.  698. 

50.  Forest,  T.W.,  and  Ward,  C.A.,  "Homogeneous  Nucleation 

of  Bubbles  in  Solutions  at  Pressures  Above  the 
Vapor  Pressure  of  the  Pure  Liquid",  J.  Chem. 

Phys .  ,  69_,  1978  ,  p.  2222. 

51.  Prausnitz ,  J.M.,  "Molecular  Thermodynamics  of  Fluid- 

Phase  Equilibria",  Prentice-Hall,  Englewood 
Cliffs,  N.J.,  1969,  p.  354. 

52.  Meyer,  C.A.£t  al . ,  eds.,  "ASME  Steam  Tables:  Thermo¬ 

dynamic  and  Transport  Properties  of  Steam" , 

ASME,  New  York,  3rd  ed . ,  1967,  p.  14. 

53.  Prausnitz,  J.M.,  op.  cit . ,  p.  148. 


54.  Frenkel,  J. ,  op.  cit.,  p.  375. 


55.  Knapp,  R.T.  et  al . ,  op.  cit. ,  pp .  55-56. 


56.  Princen,  H.M.  et  al . ,  op.  cit. ,  p.  101. 


57.  Grobner,  W. ,  and  Hofreiter,  N. ,  " Integraltaf el" , 

Springer-Verlag ,  Innsbruck,  1949,  p.  78. 

58.  Spiegel,  M.R.,  op.  cit. ,  p.  70. 


59  .  Ibid . ,  p .  Ill . 

60 .  Ibid . ,  p .  70. 

61.  Irvine,  T.F.,  Jr.,  and  Hartnett,  J.P.,  eds.,  "Steam 

and  Air  Tables  in  S.I.  Units",  Hemisphere 
Publishing,  Washington,  1976,  pp.  2-13. 

Meyer,  C.A.,  et  al.,  op.  cit . ,  p.  17. 


62. 


127 


63.  Hildebrand,  J.H.,  Prausnitz,  J.M.,  and  Scott,  R.L., 

"Regular  and  Related  Solutions",  Van  Nostrand 
Reinhold,  New  York,  1970,  p.  137. 

64.  Prausnitz,  J.M.,  op.  cit. ,  pp.  424-425. 

65.  Meyer,  C.A.  et  al . ,  op.  cit.,  p.  80q. 

66.  Defay,  R. ,  Prigogine,  I.,  and  Bellemans,  A., 

op .  cit . ,  pp.  170-171. 

67.  Himmelblau,  D.M.,  and  Arends ,  E.,  "Die  Loslichkeit 

Inerter  Gase  in  Wasser  bie  hohen  Temperaturen 
und  Drucken"  ,  Chem. -Ing . -Tech .  ,  _31,  1959, 
pp.  791-795. 


APPENDICES 


128 


APPENDIX  A 


THE  GENERAL  EXPRESSION  FOR  THE  EQUILIBRIUM  SHAPE 

OF  THE  GAS-VAPOR  BUBBLE 

The  solution  to  equation  (2.2-22),  i.e. 


dX 

dR 


-  R [ (1+e)  ~  eR2] _ 

{1  -  R2  [  ( 1+e)  -  eR2]2}''2 


(A-l) 


for  the  equilibrium  shape  of  a  gas-vapor  bubble  in  a  liquid- 
gas  solution  subject  to  a  centrifugal  field  can  be  obtained 
in  terms  of  elliptic  integrals  of  the  first  and  second  kind.^ 
Using  the  transformation 


Y 


(1  +  e) 


(A-2) 


where 


dY  =  - 2  e  Rd  R 


(A- 3 ) 


in  equation  (A-l)  results  in 


dX  = 


YdY 


r  ■»  r  3 


2e  z  [YJ  -  ( 1+e ) Y  +  e] 


(A- 4 ) 


Integrating  equation  (A-4)  gives 


X  = 


1 


7, 

2e  2 


YdY 
— n 


W 


Vr 


+  Q 


(A- 5 ) 


^The  method  of  solution  closely  follows  that  used  by 
Princen  et  al .  [56]  in  their  treatment  of  a  fluid  bubble 

immersed-!]!- and  rotating  with  a  heavier  fluid. 


129 


130 


where 

W  =  Y3  -  (1  +  e)Y2  +  e  (A-6) 

and  Q  is  a  constant  of  integration.  The  roots  Y^  >  Y^  >  Y  ^ 
of  the  cubic  term  W(Y)  are  given  by 


=  1 


j  e  +  /  e  (1  +  | ) 


1 

2  6  - 


e(l  +  |) 


For  the  interval  Y  -  Y  -  1+e  corresponding  to  0  -  R  - 
the  use  of  the  transformations 


(A-7) 


(A-  8 ) 


(A-  9 ) 

1, 


( A- 1 0 ) 

(A-ll) 


results  in  the  following  expression  for  the  integral  in 
equation  (A-5)  [57] : 

» 

YdY  _  - 2- -  [  Y  F  (K ,  ¥)  -  (Y  -Y_)  E(K,¥) 

7*  =  ^ 

4 

+  (Y1-Y3)  tanY  (1-K2  sin2?)^2] 


Y  =  Y1  -  Y2  Sin  y 

1  -  sin2'!' 


and 


K 


Y  -  Y 

_2 _ _3_ 

Y  -  Y 

1  3 


(A-12) 
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where  F (K,l)  and  E(K,1)  are  elliptic  integrals  of  the  first 
and  second  kind.  Substituting  expression  (A-12)  in  (A-5) 
gives 


X  =  1  [Y  F(K,'f)  -(Y  -  Y  )  E(K,y) 

/  e(Y1-Y3)  1  1  J 

+  (Y  -  Y  )  tanf  (1  -  K2  sin2f)  '2  ]  +  Q  (A-13) 


The  constant  of  integration  can  be  determined  by  evaluating 
equation  (A-13)  at  the  point  X=0  and  R=1 ,  where  Y=Y^. 
Equation  (A-10)  can  be  rearranged  to  give  an  expression  for 

,  i  .e . 

b 

¥  =  sin 

From  equation  (A-14)  evaluated  at  Y=Y^,  1=0.  For  1=0,  tanl , 
F  (K , 1 )  ,  and  E(K,1)  are  all  equal  to  zero.  Thus,  equation 
(A-13)  evaluated  at  X=0  and  R=1  yields 


Y  -  Y. 


1 Y  -  Y. 


(A-14) 


Q  =  0  (A-15) 

Then,  equation  (A-13)  for  the  equilibrium  shape  of  a  plane 
section  through  the  axis  of  rotation  of  the  gas-vapor  bubble 


becomes 


1 


_ _  [Y,  F(Kf¥)  -  (Yn  -  Y_)E(K^) 

/  e(Yx  -  Y3)  1  1  J 

2  9 

+  (Yx  -  Y  )  tanY  (1  -  K  sin  1)  ] 


(A- 16) 


APPENDIX  B 


DERIVATION  OF  EXPRESSIONS  FOR  THE  SURFACE  AREA,  VOLUME,  AND 

MOMENT  OF  INERTIA  OF  THE  SPHEROID 


For  the  case  of  negative  pressures  in  the  liquid- 
gas  solution,  when  e  <<  1,  the  equilibrium  shape  of  the 
gas-vapor  bubble  was  approximated  by  the  spheroid  represented 
by  equation  (3-31) ,  i.e. 


where 


(B-l) 


b  =  (1  +  e)  a 


( B-2 ) 


Expressions  for  the  surface  area.  A,  the  volume,  V,  and 
the  moment  of  inertia,  I,  of  the  spheroid  will  be  obtained 
in  terms  of  the  maximum  radius,  'a',  and  the  shape  para¬ 
meter  ,  '  e  '  . 


B . 1  Surface  Area 

If  ds  is  the  element  of  arc  length  of  the  curve 
given  in  Figure  2.2-1,  the  surface  area  can  be  expressed  as 


A 


or  using  ds 


A 


27rrds 
J  A 


4tt 


/ds i  j 

(-q— )  dX 

dx 


0 


(B. 1-1) 


(B.l-2) 
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From  equation  (B-l) 


.  2  - 1/2 
dr  _  -ax  n  x  . 

dx  “  ,2  U  ,  2] 

b  b 


(B.l-3) 


Elementary  calculus  gives  the  differential  relationship 
expressed  by 


ds 

dx 


[1  + 


(— ) 2 
dx 


V2 

] 


(B.l-4) 


Using  expression  (B.l-3)  in  (B.l-4)  results  in 


ds  _  a  n  32x2  )  /z 

dx  r  ' 1  ,  2  ; 

b 


where  the  eccentricity,  B,  is  defined  by 


2  V2 

B  =  (1  - 

b 


(B.l-5) 


(B.l-6) 


Substituting  equation  ( B .  1— 1> )  back  into  (B.l-2),  and 
rearranging  terms  yields  the  following  expression  for  A: 


A 


47TBa 

b 


dx 


(B.l-7) 


Using  integral  tables  [58]  to  perform  the  integration  in 
equation  (B.l-7)  results  in 


A  = 


4nBa 


,2 

x  ,b_ 

2  2 
Z  B 


2  2 
x  )  + 


2B‘ 


.  -1  ,Bx> 

Sln  (3} 
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(B.l-8) 


3  0 


or  upon  evaluating  the  limits, 


A  =  4Tra 


b  ,,  2.  2  b  .  -1  /Ta, 

2  (1  “  B  )  +  2b  Sln  (B) 


(B.l-9) 


The  use  of  expression  (B.l-6)  for  B  reduces  the  first  term 
inside  the  brackets  in  equation  (B.l-9)  to 


b 

2 


v2 

(1  -  B  ) 


a 

2 


(B . 1-10 ) 


For  I B I <1 ,  sin  (B)  can  be  expanded  as  a  power  series  [59], 

i .  e . 

3  5 

—  1  1  .  O  13  J 

sin  (B)  =  B  +  J  f  +  T7J  J-  +  **•  (B.l-11) 

Thus,  the  second  term  inside  the  brackets  in  equation 
(B.l-9)  becomes 


b_ 

2B 


.  -1 
sin 


(B) 


2  4 

b  bBz  3bB 

2  12  80 


(B.l-12) 


Using  equations  (B-2)  and  (B.l-6)  to  evaluate  powers  of  B, 

2 

and  neglecting  terms  of  power  e  and  greater  in  equation 
(B.l-12)  results  in 
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§3  sin  1  (B)  =  (|  +  j  e)a  (B.l-13) 

Finally,  substituting  eqautions  (B.l-10)  and  (B.l-13)  back 
into  (B.l-9)  yields  an  expression  for  the  area  in  terms 
of  ' a '  and  ' e ' ,  i . e . 


A 


4ira‘ 


(B.l-14) 


B.2  Volume 


The  volume  of  the  spheroid  can  be  expressed  as 


V  =  2 


2, 

irr  dx 


0 


(B.2-1) 


where  from  equation  (B-l) , 


2  2  . 1  x  . 

r  =  a  (1  -  — y) 

b 


(B.2-2) 


Using  expression  (B.2-2)  in  (B.2-1)  results  in 


V  =  27Ta‘ 


(1  - 


dx 

b 


0 


(B.2-3) 


Evaluating  the  integral  in  the  equation  above  gives 


T7  4  2, 

V  =  j  ira  b 


(B. 2-4) 


. 
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Substituting  equation  (B-2)  for  'b'  in  equation  (B.2-4) 
yields  an  expression  for  the  volume  in  terms  of  'a'  and  'e', 

i .  e . 

V  =  j  ira3  ( 1  +  e)  (B.2-5) 


B . 3  Moment  of  Inertia 

The  moment  of  inertia  for  the  spheroid  about  the 
axis  of  rotation  can  be  expressed  as 


(B.3-1) 


where  the  density,  pL,  is  constant.  The  integral  in  equation 
(B.3-1)  can  be  evaluated  using  cylindrical  co-ordinates  where 
'x' corresponds  to  the  cylinder  axis,  ' r'  is  the  radial 
direction,  and  0  is  the  angle  co-ordinate.  The  limits  of 
integration  for  'x' are  obtained  from  equation  (B-l)  in  terms 
of  ' r ' ,  i . e . 


x 


(B.3-2) 


For  the  volume  element  given  by 

DV  =  r  dr  d  0  dx  (B.3-3) 


expression  (B.3-1)  becomes 
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I 


fb,  2  2,  */ 

a(a  _r  > 

f  <• 

2t t 

b,  2  2.  V2 

--(a  -r  ) 

o  J 

r  dr  dO  dx 


0 


(B - 3-4) 


Performing  the  integrations  with  respect  to  'x'  and  0  in 
equation  (B.3-4)  gives 


I  = 


4TTpLb 


,2  2,  2  3  , 

(a  -r  )  r  dr 


(B.3-5) 


0 


Using  integral  tables  [60]  to  perform  the  integration  with 
respect  to  ' r'  in  equation  (B.3-5)  results  in 


I  = 


4upLb 


,  2  2, 

(a  -r  ) 

5 


-  a 


,  2  2. 
(a  -r  ) 


%  i  a 


0 


(B. 3-6) 


which  reduces  to 


I 


o  L.  4 
8ttp  ba 

15 


(B.3-7) 


Substituting  equation  (B-2)  for  'b'  in  equation  (B.3-7)  yields 
an  expression  for  the  moment  of  inertia  about  the  axis  of 
rotation  in  terms  of  'a'  and  '  e',  i.e. 


I  = 


8ttp  (1  +  e)a' 

15 
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(B.3-8) 


a 


APPENDIX  C 


NUMERICAL  DATA  USED  IN  CALCULATIONS 

C . 1  Constants 

-23 

Boltzmann's  constant  k  =  1.380  x  10  J/°K 


Avogadro ' s  number 

C.2  Specific  Volume,  vL  [61] 

Temperature 

(°C) 

0 

25 

50 

75 

100 

125 

150 

175 

200 

225 

250 

275 

300 


26 

N^  =  6.023  x  10  molecules 

kg-mole 

vL(H20) 

(m3/kg ) 

0.001000 

0.001003 

0.001012 

0.001026 

0.001044 

0.001065 

0.001091 

0.001121 

0.001157 

0.001199 

0.001251 

0.001317 

0.001404 
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C.3  Vapor  Pressure,  Pv  [62] 

The  values  of  the  vapor  pressure  used  in  the  calcula¬ 
tions  are  presented  in  Table  2  for  various  temperatures. 

These  values  were  obtained  from  the  following  equation  for 
the  vapor  pressure  as  a  function  of  the  temperature,  T,  in 
degrees  Kelvin: 


P  (0)  =  P  exp 
v  c  c 


5 

Z  (1  -  9 ) n 


9  1+K, (1-0) +K_  (1-6) 2 

D  / 


1-9 

Kg  (1-0) 2+K9 


(C . 3-1) 
T 

where  0  =  — 

c 

and  the  constants  are  as  follows: 


T  =  647.3  °K 

c 

P  =  22.120  MPa 

c 

Kx  =  -7.6912 

K2  =  -2.6080 

Kg  =  -1.6817 

K4  =  6.4233 

Kc  =  -1.1896 

Kr  =  4.1671 

D 

K?  =  2.0975 

9 

Kg  =  10 


6 
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C . 4  Partial  Molar  Volume ,  [63] 

The  partial  molar  volume  of  air  in  water  was 
approximated  by  the  value  of  nitrogen  in  water  at  a  temper¬ 
ature  of  25°C,  i.e. 


Although  rough  correlations  have  been  established 


for  the  temperature  variation  of  the  partial  molar  volume 

of  a  solute  in  a  dilute  solution  [64]  ,  in  view  of  the  low 

level  of  accuracy  involved,  it  was  decided  to  use  a  constant 

value  of  for  all  temperatures.  It  should  be  noted  that 

equation  (4-44)  for  'h'  is  not  highly  sensitive  to  the  value 

of  V? .  For  example,  at  T  =  25°C  for  PL  =  -10  MPa  and  C  =  10 
z  c 

a  fifty  per  cent  increase  in  changes  the  value  of  *h'  by 
less  than  half  of  one  per  cent. 

C.5  Surface  Tension,  y  [65] 


The  values  of  the  surface  tension  used  in  the  calcula¬ 


tions  are  presented  in  Table  2  for  various  temperatures. 

These  values  were  obtained  from  the  following  equation  for  the 
surface  tension  of  water  in  contact  with  its  vapor  as  a 
function  of  the  temperature,  T,  in  degrees  Kelvin: 


Y(T)  =  0.2358  (1  - 


1.256 


1  -  °-62511'  64735)  (N/m> 


(C.5-1) 
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The  effect  of  the  dissolved  gas  on  the  surface  tension 
was  neglected.  For  a  dilute  solution,  the  reduction  in  the 
surface  tension  of  the  pure  solvent  can  be  expressed  by 

Ayo  =  P2  RT  (C - 5-2) 

where  R  is  the  gas  constant,  and  is  the  adsorption  of 
the  solute  at  the  interface  [66] .  For  a  dilute  solution 
of  nitrogen  in  water,  the  adsorption'*'  was  estimated  to  be  of 
order  10  ^ .  Thus,  since  the  adsorption  is  very  small,  the 
effect  of  the  dissolved  gas  on  the  surface  tension  can  be 
neglected . 

C.6  Henry's  Law  Constant,  K 

The  values  of  Henry's  law  constant  used  in  the  calcu¬ 
lations  are  presented  in  Table  2  for  various  temperatures. 

For  a  dilute  solution  of  air  in  water,  can  be  expressed 

hi 

in  terms  of  K  0  and  K  _  for  dilute  solutions  of  oxygen  and 

H  A  n  J 

nitrogen  in  water. 

Consider  a  constant  temperature,  constant  pressure 
system  consisting  of  a  dilute  solution  of  air  in  water  and  an 
ideal  gas  mixture  of  air  and  water  vapor.  The  liquid  and 
gas  subsystems,  denoted  by  the  superscripts  L  and  G,  are 
separated  by  a  plane  interface.  The  argon  free  air  in  the 
gas-vapor  mixture  is  assumed  to  consist  of  twenty-one  per 

■**See  footnote  6  in  Section  2.2. 
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cent  oxygen  and  seventy-nine  per  cent  nitrogen.  The  water, 
oxygen,  and  nitrogen  are  denoted  by  the  subscripts  1,  2,  and 
3,  respectively. 

Setting  the  variation  of  the  Gibbs  potential  equal 
to  zero  for  constant  temperature,  T,  pressure,  P,  and  total 
mole  numbers,  KL  ,  yields  the  equilibrium  conditions,  i.e. 

y.L  =  y.G  i  =  1,2,3  (C.6-1) 

i  i 

Using  the  expressions  for  the  chemical  potentials  of  the 
gas  components  in  the  dilute  solution  and  ideal  gas  mixture 
in  equation  (C.6-1)  gives 

r  PX. 

y0i  (pr'T)  +  kT  ln  {~pk) 

r 

=  ip  (P ,  T)  +  kT  &  n  Ci  i=2 , 3  (C.6-2) 

Let  the  reference  system  for  each  gas  be  the  liquid  saturated 

with  the  gas  across  a  non-rigid  membrane,  permeable  only  to 

the  gas  component.  If  C.  is  the  saturation  concentration 

1o 

of  component  i,  the  equality  of  chemical  potentials  across 
the  membrane  gives 

iMP,T)  -  MoiG  (P,T)  =kTln  (^-)  i=2,3  (C.  6-3) 

i 

o 

Substituting  equation  (C.6-3)  back  into  (C.6-2)  and  choosing 
the  reference  pressure,  P  ,  to  be  equal  to  P  results  in 


1 
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PX2  = 

=  Kh2  C2  (C.6-4) 

and 


px3  . 

=  Kh3  C3  (C. 6-5) 

where  Henry's  law  constant,  K  .  ,  for  each  gas  component 

Hi 


is  given  by 

KHi  ^ 

=  i=2 , 3  (C.  6-6) 

l 

o 

Dividing  equation  (C.6-4)  by  (C.6-5)  and  rearranging  gives 


C2  = 

=  q  C3  (C.6-7) 

where 

q  = 

X2  KH3 

ST  (C.  6-8) 

X3  KH2 

For  the  air  in  the  gas-vapor  mixture,  the  ratio  of  the 
mole  fraction  of  oxygen  to  nitrogen  is  fixed,  i.e. 


X2 

X3 

0  71  'V 

=  1—1  =  0.266  (C.  6-9) 

Thus,  expression  (C.6-8) 

for  q  becomes 

q  = 

KH3 

0.266  (C.6-10) 

KH2 

If  Ku  is  Henry's 

H 

law  constant  for  a  dilute  solution 

of  air  in  water  treating  air  as  a  single  component,  then 
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P(X2  +  X3)  =  KR  (C1  +  C2)  (C . 6-11) 

Adding  equations  (C.6-4)  and  (C.6-5)  gives 

P(X2  +  X3)  =  Kr2C2  +  Kh3C3  (C . 6-12) 

From  equations  (C.6-11)  and  (C.6-12), 


KH  (C1  +  C2> 


KH2C2  +  KH3C3 


(C.6-13) 


Substituting  equations  (C.6-7)  and  (C.6-10)  into  (C.6-13) 
and  rearranging  results  in 


(1 


1.266  KH3 

KH3 

+  0.266  ) 
Ktt  o 


(C . 6-14 ) 


Equation  (C.6-14)  expresses  K  for  a  dilute  solution  of  air 

H 

in  water  in  terms  of  K  _  and  K„0  for  dilute  solutions  of 

HZ  nj 

oxygen  and  nitrogen  in  water,  respectively. 

The  values  of  KTT„  and  K  _  [67]  were  obtained  from  the 

HZ  id  J 

following  equation  for  Henry's  law  constant  as  a  function  of 
the  temperature,  T,  in  degrees  Kelvin: 


K„(Tn)  =  (1.013  x  10  )  exp 

H 


-  CCTX  +  D) 
2A 


+  D) 


4A  (BT,  +  ET. 


(MPa) 


(C . 6-15 ) 


' 
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where 


and  the  values  of  the  constants  A,  B,  C,  D,  and  E 
given  as  follows  for  oxygen  and  nitrogen: 


oxygen 

nitrogen 

A  = 

-0.0005943 

A  = 

-0.1021 

B  = 

-0.1470 

B  = 

-0.1482 

C  = 

-0.05120 

C  = 

-0.01900 

D  = 

-0.1076 

D  = 

-0.03741 

E  = 

0.8447 

E  = 

0.8510 

are 


. 
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